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SUMMARY 


This  report  covers  the  research  performed  under  AFOSR  contract  #86-0055.  The  work 
considers  the  formation  and  stability  of  partially  neutralized  plasma  clumps  that  result  when 
an  intense  electron  beam  is  passed  through  a  localized  plasma  source.  The  main  results  are: 

The  injection  of  an  intense  electron  beam  into  a  localized  gas  cloud  can  lead  to  the 
acceleration  of  a  ‘plasma  clump’  from  the  downstream  side  of  the  gas  into  vacuum,  with 
a  density  comparable  to  the  beam  density  and  at  a  speed  corresponding  to  beam  electron 
energies  and  ion  mass.  This  is  the  expected  situation  provided  the  degree  of  ionization  in 
the  gas  remains  low.  Extraction  of  the  ions  by  the  axial  electric  field  as  they  are  produced 
at  a  low  pressure  helps  maintain  this  n,  <  n*,  situation,  and  rapid  expulsion  of  newly-created 
electrons  keeps  ne  <C  n$  in  the  gas  while  there.  If  the  cathode  is  hot  enough  to 

emit  ionizing  radiation,  the  gas  should  be  shielded  (or  shadowed)  from  it  in  order  to  avoid 
photoionization. 

If,  as  at  higher  pressures,  n,-  builds  up  to  exceed  nj,,  it  is  predicted  that  electron  density 
comparable  to  nb  or  greater  can  build  up  in  the  source  gas,  leading  to  rapid  electron  heating 
by  the  two-stream  instability,  electron-impact  avalanche  of  the  remainder  of  the  gas,  and  a 
thermal  explosion  producing  a  denser  jet  of  plasma  with  much  larger  transverse  (thermal) 
energies  than  in  the  case  of  strictly  “electrostatic”  acceleration.  During  the  course  of  beam 
injection  at  pressures  on  the  order  of  1  Torr  and  above,  both  the  fast  tenuous  directed  jet 
and,  subsequently,  the  slower,  denser  jet  of  the  thermal  explosion  are  predicted,  and  are  seen 
in  a  ID  simulation.  After  beam  turnoff,  at  low  pressures,  a  plasmoid  with  n*  ~  is  seen 
in  the  simulations  to  separate  from  the  source  gas  and  propagate  downstream.  At  higher 
pressures,  simulation  shows  this  being  followed,  even  overtaken,  by  a  thermal-explosion  jet 
in  which  plasma  electrons  have  been  very  strongly  heated. 

The  interaction  of  a  propagating,  expanding  plasmoid  with  a  very  tenuous  ionized  back¬ 
ground  plasma  has  also  been  considered.  Electrostatic  effects  are  expected  to  dominate  over 
collisional  or  magnetic  effects.  Enhanced  electric  field  fluctuations  in  the  propagating  plas¬ 
moid  are  predicted  to  raise  its  transverse  electron  temperature  and  thus  contribute  to  its 
reduction  in  axially-directed  energy  flux  density  due  to  enhanced  radial  expansion. 
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1  Introduction 


This  report  considers  a  specific  system  in  which  the  formation  and  stability  of  partially 
neutralized  plasma  clumps  are  studied.  An  electron  beam  with  peak  current  of  order  20  kA, 
with  a  time  history  corresponding  approximately  to  a  25-ns  half-cycle  of  a  sinusoid,  with  a 
typical  electron  energy  of  0.5  MeV,  and  with  an  expected  final  radial  extent  of  order  1  cm, 
is  assumed  to  propagate  through  a  ~1  cm  hydrogen  gas  puff  and  create  a  virtual  cathode 
on  the  vacuum  downstream  side  of  the  gas  puff.  The  gross  beam  current  greatly  exceeds 
the  propagation  limit  for  vacuum  in  the  pipe,  so  the  beam  cannot  propagate  macroscopic 
distances  in  the  vacuum  until  it  is  largely  neutralized  there  by  ions  accelerated  forward  from 
the  ionizing  gas. 

In  what  follows,  we  describe  both  analytic  and  numerical  studies  of  the  initial  phase  of 
the  process:  the  ionization  dynamics  in  the  gas,  the  heating  of  the  resulting  plasma,  the 
formation  of  both  electrostatically  accelerated  jets  of  ions  and  thermally-produced  jets  of 
heated  plasma,  and  a  1-D  electrostatic  simulation  of  these  processes.  Then  we  describe  in 
a  preliminary  way  the  exchange  of  momentum  and  energy  between  a  propagating  partially- 
neutralized  “plasmoid”  of  finite  extent  and  low-density  background  plasma.  Some  dynamic 
simulations  of  the  “plasmoid”  formation  and  detachment  on  beam  turnoff  have  also  been 
completed,  and  we  have  made  some  analytic  conclusions  regarding  the  disturbance  of  low- 
density  background  plasma  by  the  dense  plasmoid.  Many  details  of  two-stream  interaction 
and  experimental  consequences,  and  of  electrostatic  acceleration  of  an  ion  front  from  the 
source  plasma  have  been  reported  previously  (Ref.  1  and  references  cited  therein). 

2  Source  Plasma  Ionization  -  Analytic  Theory 

In  order  to  coordinate  closely  with  our  computer  simulation  and  concurrent  experimental 
efforts,  we  have  analyzed  the  ionization  processes  in  a  ~1  cm  gas  puff  (P  <  5  torr)  through 
which  an  electron  beam  propagates  and  from  which  it  extracts  ions  to  form  a  partially 
neutralized  propagating  plasmoid  downstream.  In  describing  the  electron  and  ion  density 
buildup  by  ionization,  we  have  included  the  effects  of  electric  fields  inside  the  ionizing  gas 
puff,  and  include  ionization  by  ion-neutral  collisions.  The  gas,  puffed  transversely  across 
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the  electron  beam  path  adjacent  to  the  diode,  is  exposed  to  beam  electrons  and  cathode 
radiation  through  the  anode  aperture,  approximately  1  cm  in  diameter. 

The  simulations  realistically  model  impact  ionization  of  the  gas  by  electrons  and  ions 
but  neglect  photoionization  of  the  gas  by  cathode  radiation.  If  the  injected  beam  current  is 
much  greater  than  the  space-charge  limiting  current  for  the  system,  a  virtual  cathode  forms 
near  the  anode  plane  soon  after  the  beam  is  injected.  If  the  gas  pressure  is  high  enough,  the 
VC  moves  downstream  to  the  end  of  the  gets  region.  As  it  moves  downstream  we  find  that 
ions  are  accelerated  to  energies  several  times  the  el  beam  energy  by  coherent  motion  of  the 
ions  and  the  intense  VC  electric  fields  which  follow  the  VC.  We  also  find  that  the  potential 
in  the  gas  region  becomes  positive.  Once  the  VC  has  reached  the  end  of  the  gas  region  (if 
the  pressure  is  neither  too  high  nor  too  low)  the  system  reaches  a  quasi-steady  state  in  which 
ions  are  extracted  from  the  gas  by  the  VC  electric  fields  at  a  constant  rate;  els  are  also  lost 
from  and  propagate  downstream  at  a  constant  velocity  on  the  order  of  (eVo /mp)V2,  where 
Vo  is  the  el  beam  energy  and  mp  is  the  ion  mass.  We  also  find  that  the  secondary  els  are 
also  able  to  escape  from  the  gas  region,  even  though  the  potential  is  positive  in  that  region. 

In  the  following  sections,  we  compare  the  ionization  rates  for  el  and  ion  impact  ionization 
and  photoionization  during  the  initial  stages  of  the  beam  pulse  when  little  or  no  ionization 
is  present,  and  during  the  latter  stages  when  a  quasi-steady  state  has  been  reached. 

We  have  calculated  and  compared  (a)  impact  ionization  by  the  0.5  MeV  reflexing  electron 
beam,  (b)  photoionization  by  assumed  cathode  UV  radiation,  and  (c)  ionization  resulting 
from  collisions  of  accelerating  ions  with  neutrals.  The  beam  is  modeled  as  a  density  nj,(t) 
given  by  the  first  half-cycle  of  a  sine  wave: 

=  0.5  x  1013cm~3  sin(irt/50ns),t  <  50ns  (1) 

2.1  Electron  impact  ionization 

Assuming  the  beam  has  spread  radially  to  1  cm  diameter,  this  corresponds  to  a  peak  one¬ 
way  current  of  approximately  20  kA  for  a  radially  flat  beam  profile.  For  Hj,  we  expect  this 
primary  ionization  to  build  plasma  density  according  to 

h,(cm-3/ns)  =  .34  Ptorr  n6moI(cm-3)  sin(7r</50ns),  (2) 
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(3) 


i.e., 


n,(cm“3)  ~  5APt0rrnbmax[  1  -  cos(irt/50ns)] 
from  primary  electron  impact  ionization,  if  ion  motion  is  neglected. 


2.2  Photoionization 


The  photoionization  rate  is  essentially  independent  of  the  plasma  density  and  hence  will  be 
independent  of  gas  density  and  beam  density  during  all  stages  of  the  beam  pulse,  assuming 
that  the  cathode  radiation  is  the  same  throughout  the  beam  pulse.  We  see  that  the  pho¬ 
toionization  rate  is  strongly  dependent  on  cathode  plasma  temperature.  If  the  temperature 
is  sufficiently  high,  photoiomzation  could  produce  enough  ionization  early  in  the  beam  pulse 
to  prevent  formation  of  a  VC  in  the  gas  region.  If  the  VC  does  form  in  the  gas  region, 
photoionization  will  contribute  to  ionization  of  the  gas  and  the  VC  will  simply  move  down 
to  the  end  of  the  gas  region  sooner  and  more  quickly  than  otherwise. 

To  estimate  the  photoionization  rate,  we  first  suppose  the  cathode  plasma  to  be  a  perfect 
blackbody  radiating  as  a  several-eV  plasma.  The  radiation  flux  over  a  2x  hemisphere  solid 
angle,  from  the  15  mm  radius  cathode  (.0707  cm2),  with  unit  emissivity,  would  be 

iW)  =  1.7xl07(^)\  (4) 


with  a  more  or  less  uniform  radiation  intensity  a  distance  r  (cm)  from  the  cathode  center, 

T, 


/,(H7cm2)  *  2.7  x  106  (^) 


M 

r»2 


(5) 


The  H2  photoionization  cross-section  rises  quickly  from  the  13.6  eV  photon-energy  thresh¬ 
old  to  peak  at  8  Mbarns  =  10-17  cm2  (Ref.  2).  The  integrated  radiation  intensity 

above  14  eV  from  a  7  eV  blackbody  is  about  80%  of  the  total  flux,  so  we  could  multiply  I# 
by  0.8  and  use  the  peak  cross- sect  ion.  This  gives 


ne(cm  3)  ~ 


2.3  x  1015 
r2(cm2) 


(6) 


for  constant  cathode  conditions  during  the  first  10  ns.  (For  Tc  «  7  eV,  this  corresponds  to 
several  percent  ionization  of  the  gas,  increasing  to  ~20%  by  50  ns  when  the  beam  stops.) 
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We  refine  this  calculation  next  by  using  the  model  photoionization  cross-section  for  hy¬ 
drogen 

0  for  e*  <  6i 

<7o(ci/e^)3  for  e*  >  e, 


(7) 


tro  =  0.79  x  10-17cm2,  e,  =  13.6eV  (8) 

and  the  blackbody  spectrum  from  a  cathode  at  temperature  Tc.  From  this  we  calculate  the 
photoionization  rate  coefficient: 

‘ I—* 3  /'T  \  2 

*.(*)/>».  =  (f )  (nti/rj  (9) 

with  n  integer  and  £,(p)  the  exponential  function, 

EM  =  jf”  (10) 

The  sum  converges  rapidly  when  Te  is  not  larger  than  e,-,  and  is  shown  in  Fig.  1.  Between 
3  and  12  eV,  an  approximate  scaling  as  T 4  is  evident.  The  coefficient  multiplying  the 
sum  in  Eq.  (9)  is  approximately  2  x  109s-1.  At  cathode  temperatures  below  2  eV,  the 
sum  is  approximately  (7’c/c,)exp  —  but  the  ionization  rate  then  depends  sensitively 

on  having  the  full  blackbody  tail  of  the  cathode  photoemission  spectrum.  Given  such  a 
spectrum,  the  time  to  increase  the  ionization  fraction  fi  =  ne/na  to  1%  by  photoionization 
alone  would  be  5  ns  at  Tc  =  2.7  eV. 

Thus  we  conclude  that,  even  if  the  emissivity  of  the  cathode  is  as  low  as  0.1  and  its 
temperature  as  low  as  3.5  eV  instead  of  7  eV,  photoionization  probably  dominates  over 
impact  ionization  by  the  1  MeV  beam  electrons,  when  the  gas  is  exposed  to  hot-cathode  UV 
radiation. 


2.3  Ion-neutral  collisions 

The  modeling  of  avalanche  ionization  in  the  quasi-static  electric  fields  of  the  rising  beam 
has  been  included  in  the  analytic  theory  and  in  the  simulation.  First,  the  radial  space- 
charge  electric  field,  Er,  of  a  0.5  cm-radius  beam  with  nj,  =  0.5  x  1013  cm-3  •  sin(7r</50  ns) 
would  be  of  order  2  MV/cm  •  sin(t/16  ns)  at  the  beam  edge,  if  there  were  no  space-charge 
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neutralization.  But  in  about  10-100  ps,  photoionization  from  an  unshielded  hot  cathode 
could  produce  np  ~  tib  at  the  entrance  side  of  the  gas,  so  that  if  the  initial  plasma  electrons 
are  expelled  by  the  fields,  charge  neutralization  is  established  very  early  and  no  such  radial 
field  magnitudes  would  build  up.  Likewise,  inductive  fields  due  to  the  rise  of  net  current 
are  much  reduced  by  the  return  current  carried  by  photoionized  electrons.  But  the  residual 
electric  field  in  the  gas/plasma  may  still  be  great  enough  to  accelerate  ions  to  energies  at 
which  their  cross-section  for  ionization  of  neutrals  is  quite  large. 

To  estimate  the  ionization  due  to  accelerating  ions  of  energy  e  in  the  source  plasma,  we 
use  the  model  cross-section 

<r,(cm2)  =  lV  n  (11) 

0.73  x  10-u[e(keV)]-  798s,  e  >  43keV 

for  hydrogen.  First,  we  use  the  preliminary  estimate,  V0/£p  (diode  voltage  over  source- 
plasma  length),  for  the  average  accelerating  electric  field,  giving  e  ~  eV0z/£p  and,  to  order 
of  magnitude 

•  y  -0.3 

n,(cm'3/ns)  ~  1.63n,  Aon-  for  z/^p  >  0.1.  (12) 

This  preliminary  estimate  corresponds  to  n*  exponentiating  in  a  time 

T“(ns)  ”  °'6  [sITW^]  '  (I3) 

in  the  absence  of  ion  convection  and  other  ionization  mechanisms. 

As  ions  are  accelerated  in  the  virtual  cathode  electric  field  they  move  into  and  out  of  the 
regions  where  their  impact  ionization  cross-section  is  greatest.  To  improve  on  the  preliminary 
estimate  of  the  location  and  extent  of  hydrogen  ionization,  and  to  check  the  approximate 
validity  of  code  calculations  of  ionization,  we  made  a  simple  model  of  ion-impact  ionization 
by  ions  accelerated  from  rest  in  a  given  electric  field  distribution  that  was  chosen  to  model 
the  computed  field  distribution  in  the  gas  cloud.  A  ‘weak’  Ez  field  was  assumed  over  a 
fraction  £\f£t  of  the  gas,  and  a  stronger  Ez  field  was  assumed  over  the  virtual  cathode  layer 
of  width  £i  •<  £\  downstream  of  the  £\  region.  This  was  not  a  truly  self-consistent  calculation 
in  that  (a)  the  field  was  taken  as  given,  and  (b)  all  ion  transit-time  effects  were  ignored.  In 
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a  constant  electric  field  Et  a  proton  would  leave  the  gas  rather  quickly: 

t*,(ns)  ~  1.3[^(cm)/£2(MV/cm)]1/2  (14) 

and  on  this  timescale  we  ignore  changes  in  the  electron  beam  and  the  ionization  profile.  We 
thus  treat  the  quasi-static  problem  of  determining  the  z- dependence  of  ion  density  for  given 
electron  beam  density. 

In  this  model,  ionization  is  produced  only  by  beam  electrons  and  by  accelerated  ions. 
The  ionization  rate  at  position  z  in  the  gets  is 

♦ 

hi(z)  =  n0nb  <  aevb  >  +n0n;(z)  <  atVi  >,  (15) 

with  the  cross-section  <r,  given  as  a  function  of  energy  (or  ion  velocity)  by  Eq.  (11). 

When  the  electric  field  for  0  <  z'  <  z  has  the  constant  value  E\,  ion  energy  at  z  is  just 
e<f>(z)  —  e<j>(z')  for  ions  bom  at  z'  <  z ,  and  4>(z)  —  4>(z')  =  —  (z  —  z')E Since  the  value  of 
h,(z  =  0)  is  just  nQnb  <  aevb  >=  Sq ,  we  can  show  that  in  the  interval  z  <  ix  the  low-pressure 
equation  for  n,(z)/5o  =  R(z)  is  a  Volterra  integral  equation  of  the  second  kind, 

=  *(-§-),  (16) 

where  t  is  an  energy  variable.  Letting  u  =  €/eE\  and  w  =  z  —  u,  this  can  be  put  in  the 
standard  form 

R(z)  =  1  +  n0  /  5(z  —  w)R(w)dw  (17) 

Jo 

with  5(z)  =  (T,(e  =  xeEi). 

To  first  order  in  ~  1,  i.e.  for  low  pressures,  the  solution  of  the  integral  equation 

can  be  approximated,  using  Eq.  (11),  as 

R(z)  =  1  +  ^(noC,)xJ  +  5(n0c2)(r0-2  -  x°p2)  (18) 

for  z  >  ip,  where  for  //2, 

n0Ci  ~  516  E\ ( MV/cm) P{torr)  (19) 

n0c2  «  2.0  P(Torr)  [^(MV/cm)]-0'8  .  (20) 

and  ip  is  the  distance  required  to  accelerate  an  ion  to  the  peak  energy  of  the  ionization 
cross- sect  ion,  i.e.  43  keV  (ip  =  .215  cm  at  0.2  MV/cm).  Validity  of  the  expansion  requires 
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the  right  hand  side  of  Eq.  (18)  to  be  of  order  unity.  For  expected  Ex  fields  of  order  0.2 
MV/cm  in  the  body  of  the  gas  cloud,  the  ionization  rate  at  position  z  inboard  of  the  ‘virtual 
cathode’  at  z  =  t\  «  £,  is  given  approximately  by 

R(z)  =  1  +  2.4(P/d)  +  36(P/e,)  [(el2)a2  -  0.735]  +  0(P 2)  (21) 

times  the  electron  beam  ionization  rate,  with 

et  =  Ex/{ 0.2  MV/cm), 

0.4  <  exz  <  exlx,z  in  cm  and  P  in  Torr. 

This  result  decreases  as  the  accelerating  field  Ex  is  increased,  because  the  peak  of  the 
ion-neutral  ionization  cross-section  occurs  at  rather  low  energy,  43  keV,  compared  with  the 
energies  achieved  by  ions  accelerated  across  the  gas  by  the  E\  field  strengths  anticipated. 

One  can  add  uniform  photoionization  to  that  by  the  primary  beam  electrons  without 
difficulty,  just  by  redefining  So  to  include  both  sources.  The  quantity  R(z)  is  then  the 
ratio  of  total  ionization  to  that  by  photons  plus  beam  electrons;  i.e.,  R(z)  is  the  ionization 
enhancement  ratio  due  to  accelerated  ions. 

For  higher  pressures  (above  the  validity  of  the  first- iterated- kernel  expansion),  the  gas 
can  avalanche  due  to  ionization  by  accelerated  ions.  The  foregoing  treatment  is  limited  to 
low  pressures,  at  which  ions  leave  as  rapidly  as  they  ionize,  i.e.,  the  length  of  gas  is  less 
than  a  typical  ionization  length.  One  can  show  that  in  this  case  n<  remains  less  than  nj  if 
photoionization  is  ignored,  so  that  the  beam  is  not  charge-neutralized,  at  the  low  pressures 
for  which  Eq.  (21)  is  derived,  unless  photoionization  is  significant.  The  analysis  above  is  also 
strictly  one-dimensional  and  makes  no  provision  for  radial  loss  of  ions. 

The  mathematical  formalism  of  the  foregoing  low-pressure  model  calculation  is  insensitive 
to  the  location  or  existence  of  any  virtual  cathode  strong-field  region.  If  the  strong-field 
region  is  located  on  the  upstream  side  of  the  gas  (as  observed  in  the  computer  simulations 
at  early  times  and  low  ionization),  then  the  magnitude  of  the  term,  50,  representing  direct 
ionization  by  the  electron  beam,  is  altered,  because  beam  electrons  that  do  escape  over 
the  potential  barrier  and  stream  through  the  gets  have  lower  velocity  and  larger  ionization 
cross-section  than  would  be  the  case  with  a  virtual  cathode  at  the  downstream  end. 
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For  500  keV  -  1  MeV  electron  beams,  the  ionization  due  to  accelerating  ions  is  about  one 
order  of  magnitude  larger  than  that  by  secondary  electrons  (fief.  3),  and  it  is  the  dominant 
process  governing  the  ionization  once  photoionization  (or  any  other  process)  raises  n,-  to 
around  1013  cm-3.  However,  the  voltage  drop  across  the  source  plasma  may  collapse  when 
the  avalanche  occurs,  at  first  increasing  and  finally  reducing  the  ion  acceleration. 

As  ionization  occurs  in  the  gas,  newborn  electrons  are  at  first  rapidly  expelled  by  the 
beam  space  charge,  which  is  not  yet  fully  compensated  by  the  ions  being  formed.  This  is  the 
ion- focusing  regime  (IFR)  of  beam  propagation  through  the  gas  puff.  If  newborn  ions  are  too 
slow  to  leave  on  the  timescale  of  their  formation,  they  would  simply  accumulate  until  n,  =  n$, 
at  which  time  the  expulsion  of  newborn  electrons  would  cease.  Further  accumulation  of  ions 
by  ionization  would  then  be  accompanied  by  the  buildup  of  retained  newborn  electrons. 
This  is  important  because  it  provides  a  highly  unstable  distribution  of  electron  momenta, 
with  a  quasi-stationary  “plasma  electron”  component  in  addition  to  the  streaming  “beam 
electron”  component.  The  unstable  fluctuations  of  the  two-stream  instability  may  then  begin 
to  grow,  at  a  rate  which  depends  on  the  time- varying  densities  of  beam  and  plasma  electrons, 
and  can  reach  amplitudes  that  lead  to  considerable  heating  of  electrons,  greatly  enhanced 
electron- impact  ionization  of  the  remaining  neutrals,  and  ultimately  to  thermal  expansion 
of  the  heated  source  plasma. 

On  the  other  hand,  at  low  pressures  the  ions  created  in  the  gas  are  removed  axially  from 
the  gas  region  on  the  timescale  of  their  creation.  It  is  this  component  of  ions,  accelerated 
ahead  out  of  the  source  plasma  or  gas  into  the  beam- propagation  region,  that  allows  the 
beam  front  to  move  forward  at  the  ion  speed.  If  all  the  newborn  ions  were  extracted  at  the 
rate  they  are  produced,  the  ion  density  in  the  gas  puff  would  not  build  up  to  n,  =  n*,  and  loss 
of  IFR  behavior  would  not  occur.  The  instability  in  this  case  would  be  weak  or  nonexistent, 
and  there  is  the  possibility  that  a  plasma  object  of  beam  electrons  and  extracted  ions  might 
have  a  reasonably  low  transverse  temperature. 

For  this  reason  it  is  important  to  investigate  the  ion  creation  and  extraction  rates  to 
try  to  define  the  regions  of  parameter  space  (gas  pressure,  beam  current  density,  etc.)  over 
which  these  two  qualitatively  different  sorts  of  behavior  (IFR  and  two-stream  heating)  will 
occur. 
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One  useful  hypothesis  is  that  a  quasi-steady  state  is  reached,  in  which  the  ion  extraction 
rate  equals  the  creation  rate.  Since  all  newborn  ions  are  accelerated  in  the  quasi-static 
electric  field  of  the  beam,  if  we  know  the  time-averaged  potential  drop  over  most  of  the  gas 
region,  and  if  we  know  that  the  beam  head  is  advancing  at  a  speed 

(22) 

(with  ff,  the  beam-electron  kinetic  energy  on  entry  to  the  gas),  then  we  can  use  the  continuity 
equation  with  an  ionization  source  term  to  estimate  the  resident  ion  density  in  the  gas  puff, 
for  quasi-steady  conditions: 

drii/dz  «  Vi(z)/vi(z)  —  riiVi/v?  (23) 


with  i>i(z )  the  local  ionization  rate  (net  rate  if  recombination  is  important)  and  t),(z)  the 
electric  acceleration.  (Here  we  used  V  •  v  v/v.)  Simply  put,  if  the  volume  AbVkt  ahead  of 
the  gas  is  filled  with  ions  to  a  density  n,-  zs  n b  ( Ab  is  the  transverse  area  of  the  beam),  then 
quasi-steady  conditions  occur  when  the  spatial  integral,  over  the  gas  cloud,  of  u,(z)n0(z ),  is 
equal  to  the  ion  flux,  n^: 

<  Vi(z)n0(z)  >z  i ,  »  nbvh.  (24) 


For  a  uniform  gas  density  n0  over  a  gas  cloud  of  length  i,  then,  the  average  ionization 
rate  required  for  such  quasi-steady  conditions  is 

1/2 


(25) 


If  the  ionization  is  predominantly  done  by  ions  being  accelerated  in  the  Ez  field,  as  is  the 
case  in  the  code,  we  may  estimate,  for  a  given  Ez{z)  profile  and  quasi-static  ion  flows, 


ft. 

/  dzi/,(t;,(z,zi)), 
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with  Vi(z,z\)  =  [2 e(<j>(zi)  —  <f>(z))/rrii]l^2,<f>(z)  =  —  /0Z  Ez{z')dz\  and  z i,  the  birthplace  of  an 
ion,  as  a  parameter.  Using  the  approximate  solution  for  i/,-(z)  developed  above  for  hydrogen 
(eq.  21),  and  neglecting  ionization  created  in  the  thin  “virtual  cathode”  sheath  at  :  « 
we  have  an  approximate  condition  for  steady-state  conditions  at  low  pressure: 


5.34(F/e1){5.62(f,e1)a2-l}  =  —  -  1, 


(26) 
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(27) 


where  t\  =  f?j/(0.2MV/cm),  ia  is  in  cm,  P  is  in  Torr,  and 

k  ~  J~c  (^~)  I7,07  X  1Ql6  <Te6,(^c)l 

with  <7 eb*(0c)  the  cross-section  for  ionization  by  beam  electrons  moving  at  speed  0c.  Eq. 
(26)  is  a  condition  relating  the  gas  cloud  width  the  typical  electric  field  strength  E\  — 
(0.2MV/cm)ex,  and  the  pressure  P,  required  for  quasi-steady  conditions  (when  P  is  low 
enough).  The  electron  beam  current  has  cancelled  out  from  the  expression  because  both  the 
ionization  rate  and  the  downstream  ion  density  are  proportional  to  it,  if  photoionization  is 
absent. 

If  the  typical  E- field  value  in  the  gas  is  taken  from  computer  simulation,  then  at  low 
pressures  the  conditions  necessary  for  quasi-static  behavior  can  be  calculated  and  compared 
with  the  simulation  for  consistency. 

Ion  densities  much  higher  than  this  should  lead  to  avalanche  and  nearly  complete  ion¬ 
ization  of  the  gas.  Ion  densities  there  much  below  this  value  could  remain  quasi-static 
(extraction  balancing  production),  but  would  limit  the  density  of  the  propagating  electron 
4-  ion  beam  finger  to  the  density  available  from  ion  extraction,  i.e.,  the  full  density  nb  of  the 
electron  beam  source  would  not  be  found  in  the  propagating  plasma. 

For  beam  energies  ~  1  MeV,  k  in  Eq.  (27)  is  large  enough  that  Eq.  (26)  simplifies  to 

(28) 

for  the  hypothetical  steady  situation.  For  expected  values  of  ei,  this  would  give  values  of  P 
on  the  order  of  1  Torr.  In  fact,  the  low-pressure  approximation  is  probably  inaccurate  above 
about  0.1  Torr  (and  thus  the  denominator  expression  in  Eq.  (28)  is  not  to  be  taken  very 
seriously),  but  the  order-of-magnitude  of  P  is  probably  still  near  1  Torr,  unless  significant 
photoionization  is  present. 

3  Source  Plasma  Heating  -  Theory 

Once  the  plasma  ion  density  is  equal  to  the  electron  beam  density  (assuming  ions  are  not 
expelled  fast  enough  to  prevent  this),  the  main  collisionless  interaction  of  the  electron  beam 
with  the  ionizing  gas  puff  is  the  two-stream  instability.  This  instability  transfers  energy  from 
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the  beam  in  the  form  of  nearly  longitudinal  (electrostatic)  waves,  and  simultaneously  acts 
through  the  wave-beam  interaction  to  broaden  the  beam  in  momentum  space  as  it  advances 
through  the  plasma  region. 

If  the  electron  beam  and  background  plasma  are  relatively  cold  (i.e.,  the  momentum 
spread  of  the  beam  is  small),  then  the  energy  density  in  unstable  electrostatic  oscillations 
build  up  to  a  peak  value  of  order 

ew  =  7«»Q!  mc2^(l  +  S)~s/2,  S  <  ^  (cold  beam)  (29) 

where  a  is  a  numerical  factor  of  order  unity,  and 

S  =  /?S(n,/2n,,)''3,  (30) 


before  being  limited  by  nonlinear  processes  (Thode  and  Sudan,  ref.  7).  This  buildup  of 
wave  energy  occurs  on  a  timescale  of  several  e-foldings  of  the  instability.  The  inverse  of  the 
e-folding  time  is  approximately 

1/3 


Tx  =  Imu> 


2-y  \2 nv) 


(31) 


for  a  cold  beam  (hydrodynamic  regime  of  instability).  If  the  relativistic  electron  beam  is 
warm  or  sufficiently  scattered  in  pitch-angle  (kinetic  regime  of  instability),  the  growth  rate 
is  (Ref.  6) 


T  i  =  Imu>  « 


1 


(32) 


7nP  (A0)2’ 

where  AO  is  the  angular  width  of  the  beam  momenta. 

Ordinarily,  this  growth  is  then  limited  or  reversed  after  saturation  by  nonlinear  effects 
other  than  trapping  of  beam  electrons.  On  the  time  scale  of  ion  motion  across  a  few  Debye 
lengths,  high  electric  field  regions  (spikes)  are  formed  which  tend  to  exclude  plasma  ions  and 
electrons  by  ponderomotive  pressure  (Refs.  8-11).  This  leads  to  the  production  of  localized 
density  reductions  (cavitons).  As  these  spikes  sharpen  and  collapse  (Ref.  8),  their  Fourier 
spectrum  is  broadened.  By  this  means,  the  wave  spectrum  can  change  nonlinearly  from 
waves  with  nearly  a  single  wave  vector  to  a  spectrum  that  includes  oscillation  of  zero  phase 
velocity  and  waves  in  resonance  with  the  tails  of  the  plasma  electron  distribution.  This 
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process,  also  an  unstable  one  (Refs.  9-11)  is  referred  to  as  the  oscillating  two-stream  (OTS) 
instability. 

There  is  a  threshold  energy  for  the  onset  of  the  OTS  instability  (Ref.  9).  This  threshold  is 
reached  when  cl5  the  wave  energy  density  in  the  unstable  waves  generated  by  the  two-stream 
instability,  becomes  sufficiently  large  that  the  condition  ti/nT  ~  (&i  A#)2  is  satisfied,  where 
A £)  =  ( T/4rne 2)1^2  is  the  Debye  length  and  k\  is  the  most  unstable  wave  vector  associated 
with  the  two-stream  instability.  For  the  case  of  a  cold  beam,  Freund  et  al.  (Ref.  11)  find  the 
growth  rate  to  be 

rOTS  =  uyu/(4ir  nT)  (33) 


when  3(A:iAd)2  <  e\/nT  <  m /  mt,  and 


(34) 


when  t\ /nT  >  m/m,-  and  3(t1A/j)2. 

After  the  threshold  value  of  Ci  is  reached,  the  OTS  instability  begins  feeding  energy  to 
plasma  electrons  on  the  timescale  2/r0TS,  as  a  result  of  stochastic  acceleration  of  electrons 
crossing  the  cavitons.  Energy  is  also  fed  into  electromagnetic  waves  with  u  >  up  (Ref. 
12)  and  to  ions  (Refs.  16-18).  Because  new  beam  is  continuously  supplied  to  the  region 
throughout  the  beam  pulse,  power  is  continuously  supplied  at  a  rate  estimated  to  be 


P(e rg/cm3/s)  ~  0.4  x  10 


is  ( _ ?L 


\1012cm" 


(35) 


Some  significant  fraction  (~  1/2)  of  this  power  goes  into  heating  of  plasma  electrons,  which 
subsequently  lose  energy  by  inelastic  collisions,  by  expansion  (Refs.  19-22),  or  by  loss  from 
the  system.  The  energy  input  rate  per  plasma  electron  corresponds  to  a  “heating  rate” 

Te  -I-  — (T,  +  2e,/3)  »  170/(1  4-  S)-5/2(l03nb/ri')5/3  eV/ns,  (36) 

ne 

where  /  is  the  fraction  of  deposited  power  delivered  to  electrons,  and  e,  is  the  ionization 
energy  per  atom.  Cooling  by  electron  loss  at  their  thermal  speed,  by  comparison,  occurs  at 
a  rate 

re(eV/ns)  -  1330[re(keV)]3/%-\  (37) 
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for  a  plasma  of  effective  size  lv  (cm).  And  classical  electron  heat  conduction  could  lead  to  a 
cooling  rate 

T,(eV/ ns)  ~  -380[re(keV)]7/2^2,  (38) 

with  an  effective  temperature  gradient  scalelength  It  (cm),  if  the  diode  or  surrounding  gas 
serves  as  a  heat  sink. 

For  1-centimeter  scalelengths,  then,  one  might  expect  mean  electron  thermal  energies 
roughly  of  the  order  of  seven  hundred  eV  in  the  source  plasma,  with  possible  nonthermal 
‘tails’  on  the  electron  energy  distribution  there  extending  out  to  many  keV.  Once  n,-  exceeds 
nj,  a  10  ns  beam-rise  represents  a  few  thousand  plasma  oscillation  periods  (their  variation 
included),  or  several  hundred  growth  times,  and  one  may  thus  expect  a  fully-developed  two- 
stream  instability  in  the  source  plasma  unless  gradients,  quasi-static  electric  fields,  or  other 
effects  significantly  damp  the  instability.  (The  plasma  is  only  about  20  wavelengths  thick, 
at  the  most  unstable  wavelength.) 

Such  a  microinstability  probably  has  three  important  effects  besides  generation  of  mi¬ 
crowave  diagnostic  radiation:  it  can  significantly  heat  the  plasma,  leading  to  rapid  electron 
avalanche,  it  can  make  the  plasma  electrons  effectively  more  collisional  and  less  tied  to  mag¬ 
netic  field  lines  in  laboratory  experiments,  and  it  can  spread  the  energy  distribution  of  the 
beam  electrons.  The  first  and  last  are,  in  fact,  observed  in  the  simulations  described  in  the 
next  section. 

Some  stabilizing  influences  which  could  reduce  the  plasma  heating  due  to  this  microin- 
stability  have  been  studied.  They  are: 

(a)  electromagnetic  effects,  which  can  lead  to  a  threshold  current  density  for  onset  of  the 
instability  growth  (Refs.  23-25); 

(b)  radial  and  longitudinal  boundary  conditions  (Refs.  24-27); 

(c)  axial  density  gradient  in  the  source  plasma  (Refs.  2S-35); 

(d)  virtual-cathode-related  electric  field  effects  on  instability  growth  (Refs.  36,  37); 

(e)  electron  reflexing  from  a  more  realistic  (non  1-D)  beam  head. 
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3.1  Threshold  currents 


These  thresholds  (Ref.  23)  were  found  to  be  much  lower  than  the  currents  of  our  simulations 
and  the  University  of  Maryland  experiments.  Instability  growth  should  thus  be  unaffected. 

3.2  Radial  and  longitudinal  boundary  conditions 

The  effects  of  radial  boundaries  were  considered  in  Refs.  (24)  and  (25).  In  anticipated 
experiments,  the  beam  and  pipe  were  large  enough  to  contain  several  wavelengths  of  the 
most  unstable  wave.  Stabilization  by  finite-radius  effects  is  thus  not  expected. 

Longitudinal  boundary  effects,  including  wave  reflection,  can  have  much  more  impact 
(Refs.  26,  27). 


3.3  Axial  finite-length  effects  and  ion-density-gradient  effects  on  two-stream 
growth 


Using  calculations  of  the  spatial  exponentiation  length  of  the  two-stream  instability  published 
by  Rowland  (Ref.  26),  we  find  the  growth  length  in  a  uniform-density  plasma  with  >  nh 
to  be 

^  (cm)  ~  0.5  x  10-2Te'/.3(7  -  l)1'6  ^10-^C-  ^  '  .  (39) 

Thus,  even  as  the  plasma  electrons  become  heated  to  ~10  keV,  a  1  cm- long  plasma  source 
region  always  provides  20  or  more  growth  lengths  of  the  instability  at  np  =  1015  cm-3. 
nb  =  1013  cm-3,  7  =  2,  provided  np  is  uniform. 

Next,  from  the  work  of  Schmidt  (Ref.  28),  we  have  an  estimate  of  the  maximum  number 
of  e-foldings  (of  the  unstable  fields)  allowed  by  any  given  density  gradient  in  n^: 


AT  =  0.8  x  10’A  (— -V  (40) 

\1015  cm-3/  V10l2cm-3/ 

(for  the  non-collisional  cold-beam  instability),  where  h  is  the  density  gradient  scaleheight, 

h~X  =n~Z\dnP'ldz\.  (41) 


2/3 


The  weak  dependence  of  Eq.  (38)  on  7  has  not  been  included  here  because  the  theory  is 
nonrelativistic.  For  anticipated  parameters  as  above,  then,  nonlinear  effects,  rather  than 
density  gradients,  would  limit  the  e-folding  if  the  beam  is  cold. 
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3.4  Effect  of  axial  electric  field  on  instability  growth 

When  the  residual  quasi-static  field  in  the  plasma  is  of  the  order  of  the  beam-electron  trap¬ 
ping  field  Et  and  is  in  the  direction  to  retard  beam  electrons,  a  significant  damping  of  the 
instability  can  occur  (Refs.  36,  37). 

3.5  Qualitative  effect  of  radial  beam  electron  loss  instead  of  axial  reflexing 

In  the  code  of  section  4,  which  describes  purely  axial  motion  of  beam  and  plasma  electrons 
as  if  they  were  fully  magnetized  and  did  not  transfer  energy  to  perpendicular  motion  or 
temperature,  beam  electrons  reaching  the  “beam  head”,  i.e.  the  advancing  ion  front,  reflex 
axially  back  toward  the  source  (except  for  a  few,  which  may  leak  forward  and  contribute  the 
limiting  current  propagating  in  the  vacuum).  However,  in  a  not-strongly- magnetized  real  sit¬ 
uation,  most  of  those  beam  electrons  reaching  the  front  would  be  lost  radially  -  some  sooner, 
some  later  -  rather  than  specularly  reflected  from  the  front  axially  back  to  the  source.  Thus, 
code  physics  dependent  on  the  presence  of  the  counterstreaming  reflected  beam  electrons  is 
probably  not  to  be  taken  as  representative  of  real  experimental  conditions.  Such  physics  in¬ 
cludes  heating  of  the  propagating  plasma/beam  by  two-stream- unstable  waves  arising  from 
the  beam-electron  counterstreaming.  It  also  includes  a  factor  of  two  in  the  propagating 
region  electron  density  and  in  the  beam-electron  density  in  the  gas/source  region. 

Radial  loss  of  beam  electrons  at  the  head  of  the  propagating  plasma  also  may  alter  the 
depth  of  the  beam  potential  well  there,  because  of  the  negative  space  charge. 

4  Propagating  Plasma  Jet  Formation  -  Theory 

It  appears  that  ions  can  be  accelerated  from  the  plasma  by  the  motion  of  the  steep  potential 
gradient  at  the  virtual  cathode  and,  when  applicable,  by  the  unstable  electrostatic  waves 
brought  to  large  amplitudes  by  streaming  instability  of  the  beam  with  the  source  plasma 
electrons  (Refs.  16-18).  These  accelerated  ions  propagating  into  the  low  density  background, 
together  with  the  (reflexing)  electrons  of  the  beam,  can  form  a  low-density,  partially  neu¬ 
tralized  plasmoid,  which  has  an  axial  extent  determined  by  the  beam  pulse  duration,  the 


15 


acceleration  time,  and  the  accelerated  ion  velocity, 

~  (Tei  ~  Ta)Vi.  (42) 

The  density  of  this  plasmoid  is  <  neb,  the  accumulated  density  of  reflexing  beam  electrons. 

Another  mechanism,  qualitatively  different,  also  produces  a  propagating  dense  plasma 
when  the  source  plasma  is  sufficiently  long  in  the  axial  direction:  the  unstable  waves  not 
only  accelerate  a  few  ions  to  quite  high  velocities  (t>/,  ~  /9c)  but  also  can  heat  a  much 
larger  fraction  of  source  plasma  electrons,  which  causes  the  source  plasma  to  expand  (the 
hot  electrons  electrostatically  drag  the  plasma  ions),  and  this  thermal  expansion,  somewhat 
contained  in  other  directions,  forms  a  plasma  jet  in  the  beam  direction  which  can  have 
density  n„  »  neb  but  has  expansion  front  speed  vf  ~  Me,  <  /3c  (M  =  Mach  number,  c,  = 
sound  speed  based  on  heated  source-plasma  electron  temperature  —  “tail”  temperature  if 
appropriate.) 

This  plasma  jet  also  forms  a  propagating  plasmoid,  terminating  its  ejection  from  the 
source  plasma  when  the  latter  either  cools  or  is  exhausted.  The  cooling  and  evacuation  of 
the  source  plasma  may  both  be  related  to  the  instability  heating  in  a  complex  nonlinear 
way:  the  instability  growth  should  be  altered  by  source-plasma  density  gradients,  and  the 
plasma  heat  source  is  the  thermalizing  of  energy  from  the  unstable  waves.  For  the  warm 
beam  case  (which  our  one-dimensional  computer  simulations  indicate  applies  after  a  few 
ns)  the  coupled  instability  and  hydrodynamics  was  treated  in  Refs.  19-22;  the  results  are 
dependent  on  the  plasma  thermal  conductivity  and  on  the  sensitivity  of  the  Langmuir  wave 
growth  rate  to  the  density  gradients  resulting  from  expansion  (Refs.  23-31).  The  instability 
growth  is  also  affected  by  the  steep  quasi-static  potential  gradient  (Refs.  32,33)  so  that  the 
two  mechanisms  are  not  completely  independent. 

Thus  after  a  time  lag  it  is  possible  that  two  plasma  fronts  will  be  seen  advancing  down¬ 
stream:  one  for  a  fast,  tenuous  plasma  ejected  electrodynamically  and  one  for  a  slower  denser 
plasma  ejected  thermally.  The  timescale  of  our  computer  simulations  is  long  enough  to  see 
the  formation,  if  not  the  termination,  of  the  faster  one,  but  also  shows  the  beginning  of  the 
thermally-produced  one  when  the  temperature  reaches  ~  1  keV.  The  energies  in  the  fast  and 
slow  plasmoids  can  be  estimated: 
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Fast:  n  ~  r»ej 


^ £pnmiV 2  ~  ^(re4  -  ra)ne4ro,t;3,  (43) 

Slow:  n  <./*  -  n,*,/(c.rc) 

~  ^^ne4(ne4/nc)1/3(72  -  l)mec2,  (44) 

i  b 

in  ID,  where  vf  ~  (me/mi)^2<P,  neb  is  the  reflexing  electron  beam  density,  and  it  was 
assumed  that  the  heated  source  plasma,  with  “initial”  density  nt  and  length  t„  receives  only 
a  fraction  (ne4/n,)^3( 7+1)  of  the  beam  kinetic  energy  density,  1/3  of  this  going  into  directed 
kinetic  energy  via  thermal  heating,  and  undergoes  a  free  expansion  (here  one-dimensional) 
on  timescale  ra. 

The  ~  (rather  than  ~)  refers  to  the  fact  that  once  the  reflexing  electron  beam  is  “warm”, 
the  growth  rate  of  unstable  energy  deposition  drops  and  the  saturation  energy  density  is  also 
reduced,  so  that  less  energy  may  be  available  as  thermal  to  drive  the  ablation  of  the  slow 
plasma  jet. 

Numerically  it  is  reasonable  to  expect  the  fast  low-density  plasmoid  to  contain  more 
energy  than  the  slower  denser  one  if  n  ~  ne4  <  n,  and  v,  ~  0.1  c. 


"Wcm  J>~14(lo^)(l^)(o^)3  <45) 

Slow: 

"('/cm”)  ~  ™  ><  ><>-3  (l o^r  «P^1'3  fe)  •  (46» 

with  fi  the  degree  of  ionization  in  the  source  plasma  and  the  length  over  with  beam 
energy  is  (partly)  deposited. 

However,  if  the  electron  beam  is  made  to  deposit  most  of  its  energy  in  the  source  plasma 
rather  than  leaving  only  an  energy  density 

E2/8it  ~  (7  +  l)(ne4/n,)1/3(7  -  1  )nef>mec2,  (47) 

or  if  this  wave  energy  density  can  be  transformed  rapidly  enough  into  thermal  energy  by 
nonlinear  or  collisional  process,  an  “anomalous  intense  driver”  (Ref.  38)  of  a  more  energetic 
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denser  plasma  jet  may  develop  with 

CV/cm2)  ~  5(7  -  Dntmjffn  ~  100<7  -  1)  (10,£m.,)  (j^)  •  <48> 

(This  may  require  a  longer  plasma.) 

The  relative  merits  of  such  a  hot  plasma  jet,  as  compared  with  the  tenuous  but  energetic 
electron  and  accelerated  ion  plasmoid,  should  now  be  considered.  Specifically,  we  consider 
transverse  temperature.  Collectively  accelerated  ions  can,  in  principle,  have  very  little  trans¬ 
verse  temperature,  but  a  plasma  jet,  even  after  expansion  cooling,  would  be  expected  to  have 
transverse  beam-heated  thermal  energy  comparable  with  its  directed  energy  (Ref.  39).  In  the 
absence  of  an  axial  confining  magnetic  field,  or  of  transverse  cooling  by  some  rapid  means, 
radial  expansion  velocities  would  be  only  somewhat  less  than  axial  expansion  velocity.  For 
a  1  keV  hydrogen  plasma  this  would  be  of  order  108  cm/s.  Initially,  at  least,  such  a  dense 
plasmoid  would  be  quasi-neutral  and  have  J,  %  \Je\  and  so  would  not  be  constrained  by 
magnetic  self-focusing  due  to  its  own  current.  By  contrast,  the  lower  density  plasmoid  with 
accelerated  ions  need  not  be  quasi-neutral,  but  would  have  electrostatic  forces  from  its  ini¬ 
tial  space  charge  imbalance  and/or  from  the  preferential  scattering  of  one  species  (usually 
electrons)  out  of  the  co-moving  plasma  by  collisions. 

The  radial  forces  on  electrons  and  ions  in  a  partially-neutralized  long  cylindrical  beam 
(plasmoid)  of  “co-moving”  electrons  and  ions  can  best  be  considered  in  the  plasmoid  frame, 
once  this  object  is  detached  from  its  source  and  propagating  in  a  tenuous  background.  This 
analysis,  although  the  electrons  may  have  relativistic  “temperature”1  in  this  frame,  transforms 
back  to  the  laboratory  frame  rather  simple  because  the  relative  velocity  is  less  than  about 
0.1  c.  In  the  plasmoid  frame  there  is  almost  no  axial  net  current  (except  possibly  for  fast 
oscillations  with  zero  average,  which  should  damp  rapidly  by  phase  mixing).  For  times  short 
compared  with  electron  isotropization,  the  transverse  electron  temperature  could  be  assumed 
nonrelativistic  as  well,  but  the  rapid  growth  of  high-frequency  fluctuations  in  the  formation 
state  (many  with  k±  ^  0)  indicates  that  a  better  assumption  might  be  nearly- isotropic 
electrons  with  most  of  their  initial  directed  energy  now  appearing  as  thermal  energy  in  the 
plasmoid  frame  (Ref.  39).  Ignoring  the  background  ions,  then,  this  long  plasmoid  expands 
radially  under  its  pressure  gradient  and  net  space-charge  forces.  While  at  early  times  the 


electron  radial  accelerations  may  significantly  exceed  those  of  the  ions,  in  a  few  plasma 
periods  the  radial  space-charge  oscillations  damp  to  an  ambipolar  one-fluid  expansion  at  a 
few  times  the  ion  sound  speed.  Provided  the  electrons  do  not  escape  at  early  times  enough 
to  lower  Wpe  below  u we  may  take  the  asymptotic  time  dependence  of  u>pe(t)  as  controlling 
the  damping  of  the  radial  space-charge  pulsations. 

For  radial  expansion  at  the  sound  speed  ~  (Te/m,)1/2,  the  ratio  of  radial  to  axial  speeds 
in  the  laboratory  frame  is  roughly 

t  ~  ((7  -  A7  -  l)m/m,j1/2^2_1(ro/r)2,  (49) 

where  A7  is  the  energy  loss  in  the  formation  and  counterstreaming  stages  of  early  develop¬ 
ment,  in  units  of  me2;  m  without  subscript  refers  to  the  electron  rest  mass,  and  adiabatic 
cooling  (Te  oc  r~2)  has  been  assumed.  For  the  experimental  beams  in  our  laboratory  (and 
simulations),  with  7  as  2  and  A7  as  0.1,  mt  =  1836  m,  and  0Z  ~  0.02,  this  would  give  an 
initial  expansion  cone  with  initial  half-angle  of  about  49°  or  larger,  which  is  unacceptable 
for  some  applications,  and  is  probably  made  worse  by  increasing  7. 

Ignoring  electrostatic  assistance  of  the  expansion,  one  derives 


r(z)  ~  ro[l+3KZ0z  Vq1]^3, 


(50) 


with 

k  ~  [(7  -  A7  -  l)m/m,]1/2.  (51) 

For  the  laboratory  conditions  above,  this  would  lead  to  expansion  from  r0  =  0.5  cm  to 
r  ~  9  cm  in  1  m  of  travel  at  0Z  =  .02  (for  protons). 

The  conclusion  of  this  portion  of  the  analysis,  then,  is  that  only  by  using  higher-mass 
ions  and/or  preventing  electron  isotropization  can  the  radial  expansion  rates  in  vacuum  be 
kept  low  enough.  It  is  unlikely  that  electron  isotropization  can  be  prevented,  as  it  results 
from  the  electron  reflection  at  the  advancing  plasmoid  front  as  well  as  from  unstable  waves 
with  radial  components  of  k. 

The  low-density  background  plasma  streaming  through  this  plasmoid  in  its  frame  can 
affect  the  radial  expansion  if  its  electrons  are  slowed  more  than  its  ions  (generating  a  net 
current  and  thus  an  aximuthal  magnetic  field  in  the  plasmoid  frame);  the  magnetic  effect  is 
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slight  if  the  magnetic  energy  density  is  small  compared  with  the  thermal  energy  density  of 
the  isotropized  plasmoid  electrons.  This  is  generally  the  case  until  the  piasmoid  radius  is 
quite  large 

iL  <  |eV[A (nb/3)]2  (52) 

where  A(nb/f?)  is  the  value  of  —  nbe0e  in  the  background  (in  the  plasmoid  frame)  caused 
by  effective  “collisions”  with  the  denser  plasmoid,  and  r  is  the  radial  extent  of  the  plasmoid. 
Thus, 


nvTp 

52/8* 


2 np  (7-  A7-  ljmc2^ 
*r2e2  [A(nb/?)]2  r2 


(53) 


Here  c/up  is  the  collisionless  plasmoid  skin  depth,  and  we  have  allowed  the  plasmoid  to  cool 
by  adiabatic  expansion  from  its  creation  radius  r0.  Taking  the  plasmoid  density  reduction 
(from  radial  expansion)  into  account,  the  magnetic  effects  are  important  only  once 


~  2.8  (— ) 
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h  -  A7  - 


(54) 


where  A(nb^)/np  is  approximately  constant,  and  where  the  skin  depth  at  injection  is 

Il3\  J/2 


c 

ujpo 


-«(5) 


(cm). 


(55) 


Since  np/nb  is  quite  large  (nb  here  is  the  background  electron  density  in  the  “vacuum”  region), 
magnetic  braking  of  the  plasmoid  expansion  should  occur  only  for  large  plasmoid  radius,  if 
then. 


If  background  electrons  are  scattered  into  a  halo  surrounding  the  plasmoid  (leaving  the 
background  ions  in  the  plasmoid),  there  is  a  slight  electrostatic  inward  tracking  force  which 
partially  offsets  any  electrostatic  outward  force  on  plasmoid  electrons  due  to  their  own  charge 
excess  (i.e.,  nep  —  n<p  >0). 


5  ID  Simulation  of  Source  Plasma  and  Plasma  Jet  Formation 

[Basic  PIC  code  and  results  are  published,  J.  Appl.  Phys.  67,  1650  (1990).  A  copy  is 
attached.] 
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5.1  Beam  Cutoff  Studies 

To  investigate  the  possibility  of  forming  a  charge-  and  current-neutral  plasmoid  in  our  system, 
we  modified  the  PIC  code  which  we  have  been  using  to  simulate  plasma-jet  formation  so 
that  the  electron  beam  cuts  off  in  the  middle  of  a  run.  Results  are  shown  for  one  such 
run  in  which  a  500-kV,  15-kA,  5-mm- radius  electron  beam  is  injected  into  a  4-cm-radius, 
40-cm-long  drift  tube  with  a  600-mTorr,  i-cm-wide  cloud  of  hydrogen  gas  located  next  to 
the  anode  plane. 

In  Figure  1,  the  density  profile  for  each  particle  species  and  the  net  charge  density  profile 
are  shown  at  intervals  of  2000  time  steps  (which  corresponds  to  about  1.5  ns  of  real  time) 
after  the  beam  has  been  turned  off.  Before  the  beam  is  turned  off,  a  plasma  jet,  composed 
of  beam  electrons  and  plasma  ions  with  nearly  equal  densities  of  approximately  1012  cm'3 
extends  from  the  end  of  the  source  plasma  to  a  distance  approximately  20  cm  downstream. 
(The  source  plasma  can  be  identified  by  the  large  spike  in  the  plasma  ion  density  profile.) 
After  the  beam  is  turned  off,  the  tail  of  the  plasma  jet  detaches  from  the  source  plasma, 
forming  a  plasmoid  composed  of  plasma  ions  and  beam  electrons.  After  detaching  from  the 
source  plasma,  the  plasmoid  continues  downstream  at  a  velocity  of  between  0.04c  and  0.05c. 

As  time  progresses,  the  number  of  particles  in  the  source  plasma  continues  to  increase, 
even  though  the  electron  beam  has  been  turned  off.  At  the  same  time,  L.*e  source  plasma 
begins  to  expand  downstream.  Unfortunately,  the  combination  of  these  two  effects  makes 
simulation  of  the  plasmoid  difficult  at  later  times.  Because  of  computer  limitations,  the  code 
halves  the  number  of  macroparticles  in  the  simulation  whenever  the  number  becomes  too 
large.  As  the  number  of  particles  in  the  source  plasma  increases,  the  number  of  macropar¬ 
ticles  in  the  plasmoid  is  therefore  constantly  being  reduced.  Also,  as  the  source  plasma 
expands  outward,  it  overtakes  the  plasmoid  and  the  plasmoid  is  no  longer  distinguishable 
from  the  source  plasma. 

Figure  2  shows  the  density  profile  for  each  particle  species  for  a  run  in  which  the  source 
plasma  was  artificially  removed  after  the  plasmoid  had  detached  itself  from  the  source  plasma. 
The  density  profiles  are  shown  20,000  time  steps  (15  ns)  after  the  beam  was  turned  off. 
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5.2  Pressure  Variation  Studies 

We  have  performed  a  number  of  simulations  with  different  neutral  gas  pressures.  Quali¬ 
tatively,  we  were  able  to  identify  four  pressure  regimes  -  ultra- low.  low,  intermediate,  and 
high. 

In  the  ultra-low-pressure  regime,  the  neutral-gas  pressure  is  too  low  to  produce  a  signif¬ 
icant  amount  of  ionization  during  the  simulation.  No  significant  movement  of  the  virtual 
cathode  is  observed,  the  electron  current  which  is  measured  downstream  equals  the  value  of 
the  space-charge  limiting  current  for  the  system,  and  no  ions  are  detected  downstream. 

In  the  low-pressure  regime,  sufficient  ionization  is  produced  to  cause  motion  of  the  virtual 
cathode  downstream  through  the  gas  cloud.  As  the  virtual  cathode  moves  through  the  gas, 
ions  from  the  ionized  gas  are  accelerated  to  energies  several  times  the  electron  beam  energy. 
The  potential  in  the  region  behind  the  virtual  <  uhode  also  becomes  positive.  The  gas 
pressure  is  too  low,  however,  to  maintain  the  forward  position  of  the  virtual  cathode,  and 
the  virtual  cathode  quickly  returns  to  its  original  position  near  the  anode  plane.  The  electron 
current  measured  downstream  slightly  exceeds  the  space-charge  limiting  value,  and  a  weak 
ion  current  is  detected  downstream. 

In  the  intermediate- pressure  regime,  the  gas  pressure  is  sufficiently  high  that  the  virtual 
cathode  is  able  to  maintain  its  position  after  reaching  the  end  of  the  gas  cloud,  and  a 
quasi-steady  state  is  achieved.  In  the  steady  state,  ions  are  created  at  a  rate  which  just 
balances  the  rate  at  which  they  are  extracted  from  the  source  plasma  by  the  electric  fields 
of  the  virtual  cathode,  which  is  located  just  downstream  of  the  source  plasma,  and  the  rate 
at  which  they  are  ejected  from  the  source  plasma  through  the  anode  plane  by  the  positive 
potential  in  the  source  plasma.  The  plasma  electrons  are  also  ejected  from  the  source  plasma 
through  the  anode  plane  by  a  mechanism  which  we  have  not  yet  identified.  The  majo  ity  of 
beam  electrons  continue  to  be  reflected  back  through  the  anode  plane  by  the  virtual  cathode. 
However,  a  significant  fraction  travel  past  the  virtual  cathode  and  down  the  drift  tube  until 
they  reach  the  head  of  the  beam  of  ions  which  have  been  extracted  from  the  plasma  by  the 
virtual  cathode  fields.  Most  of  the  electrons  are  reflected  at  the  head  of  the  beam  of  ions  but 
some  continue  downstream,  creating  a  current  equal  to  the  space-charge  limiting  current. 

In  Figure  3,  we  show  the  electron  and  ion  momentum  distributions  at  a  point  in  the 
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plasma  jet  downstream  from  the  source  plasma  but  upstream  from  the  head  of  the  jet  for 
one  run  in  the  intermediate  pressure  regime. 

In  the  high-pressure  regime,  the  gas  pressure  is  too  high  to  maintain  a  quasi-steady  state. 
The  number  of  ions  and  electrons  in  the  source  plasma  increases  with  time,  and  the  electron 
temperature  in  the  source  plasma  also  increases  with  time.  At  a  certain  time,  the  plasma 
electrons  become  sufficiently  energetic  that  they  suddenly  jet  downstream,  dragging  the  ions 
in  the  source  plasma  downstream  with  them. 

6  Momentum  and  Energy  Exchange  Between  a  Dense  Propagat¬ 
ing  “Plasmoid”  and  a  Tenuous  Background  Plasma 

We  have  analyzed  here  to  a  limited  extent  both  collisional  and  electrodynamic  processes. 
Under  almost  all  circumstances  the  latter  dominate  when  the  background  is  tenuous.  The 
collision  processes  examined  included  large-angle  scattering  and  charge  exchange  (for  weakly- 
ionized  hydrogen  background)  between  plasmoid  particles  and  background  particles.  More 
significant  are  the  effects  on  the  energy  and  momentum  exchange  mediated  by  a  level  of 
fluctuations  enhanced  above  the  thermal  level  and  mediated  by  the  non-neutrality  of  the 
plasmoid  when  its  net  charge  density  exceeds  the  background  plasma  density. 

Electron  scattering  within  the  dense  plasmoid,  in  its  own  frame,  continually  repopulates 
the  electron  Maxwell  tail  above  the  electrostatic  potential-well  depth,  allowing  electrons  to 
escape.  The  potential  well  depth  adjusts  so  that  the  ion  expulsion  rate  due  to  the  potential 
equals  the  electron  loss  rate.  When  Coulomb  collisions  dominate  the  scattering  loss,  this 
ambipolar  potential,  neglecting  the  background,  is  e^/Te  ~  occurring  over  the 

density  gradient  scalelength  of  the  dense  plasma  (roughly  half  its  dimension). 

The  effects  of  the  low  density  streaming  background  are  (if  4>  >  0): 

(a)  Background  electrons  ahead  of  the  plasmoid  are  radially  attracted,  weakly  magnetized 
by  their  net  current  (~  enb Au„),  and  axially  slowed  temporarily,  thus  increasing  their 
density  locally.  If  npi  —  n ^  is  larger  than  the  resulting  neb  —  n.j,  this  is  of  little  consequence 
for  plasmoid  electrons. 

(b)  In  the  frame  of  the  dense  plasma,  the  flowing  low-density  background  plasma  contains 
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a  (quasi-static?)  nonneutral  downstream  wake.  While  background  ions  Me  too  heavy  to  have 
their  motion  strongly  perturbed,  for  a  positively-charged  plasmoid  the  background  electrons 
are  accelerated  radially  into  and  out  of  the  plasmoid,  being  temporarily  ‘heated’  in  the 
process.  These  accelerations  can  occur  on  the  electron  inertial  timescale,  i.e.,  u?^1,  thus 
establishing  a  plasma  space-charge  oscillation  (or  standing-wave  perturbation)  which  is  the 
large-amplitude  version  of  the  perturbation  of  background  plasma  by  a  test  charge.  Here  the 
‘test  charge’  has  finite  extent  and  ecf>  ~  7^,  so  the  space-charge  perturbation  in  the  wake 
can  have  Sne/nei  ~  1  and  e<f>  ~  7)*  over  a  radial  scalelength  comparable  to  the  dense-plasma 
radius. 

(c)  By  enhancing  the  electric  field  fluctuation  level  in  the  plasmoid,  the  streaming  back¬ 
ground  enhances  the  effective  energy-exchange  collision  rate  of  the  dense  electrons  above  the 
Coulomb  value  t/e;  thus  the  energy  tail  of  the  electron  distribution  function  would  be  popu¬ 
lated  faster,  assuming  0c  >  2 vth  in  the  plasmoid.  In  order  to  compensate  for  the  increased 
escape  of  electrons  over  the  electrostatic  barrier,  the  plasma  potential  must  be  increased  to 
allow  both  ions  and  electrons  to  be  lost  at  equal  rates.  The  potential  increases  by  an  additive 
(.n(vlfvi),  where  vL  is  the  electron  loss  rate  in  the  absence  of  enhanced  fluctuations  and 
v'L  is  the  enhanced  loss  rate.  In  steady  state,  vl  is  proportional  to  i/(.  Once  u*  vt  the 
potential  change  should  scale  with  the  logarithm  of  the  fluctuation  level,  and  the  energy  in¬ 
vested  in  quasi-static  potential  is  proportionately  larger.  The  fields  are  thus  larger  for  given 
plasmoid  size,  when  0c  vth,  but  the  energy  dissipation  rates  are  not  strongly  affected  by 
these  electrostatic  fields. 

However,  the  fluctuations,  supplied  by  background  streaming  energy  (i.e.,  in  the  lab 
frame  by  plasmoid  slowing),  also  serve  to  heat  the  plasmoid  electrons  until  0c  ~  Then 
the  loss  region  (e  >  e<f>)  is  depopulated  by  trapping  background  electrons  almost  as  fast  as 
it  is  populated  by  tail-filling  from  the  denser  electrons  (as  one  can  see  by  looking  at  phase- 
space  orbits  near  the  trapped-untrapped  boundary),  and  the  plasma  potential  and  associated 
quasi-static  E- field  decrease. 

(d)  The  magnetic  field  in  the  plasmoid  frame  due  to  the  change  in  vz  of  the  streaming 
background  electrons  slightly  impedes  the  thermal  expansion  of  the  plasmoid  in  the  radial 
direction,  as  discussed  above. 


24 


7  Discussion  and  Conclusions 


The  injection  of  an  intense  electron  beam  into  a  localized  gas  cloud  can  lead  to  the  accelera¬ 
tion  of  a  ‘plasma  clump’,  with  density  comparable  to  the  beam  density,  from  the  downstream 
side  of  the  gas  into  vacuum,  at  a  speed  corresponding  to  beam  electron  energies  and  ion  mass. 
This  is  the  expected  situation  provided  the  degree  of  ionization  in  the  gas  remains  low.  Ex¬ 
traction  of  the  ions  by  the  axial  electric  field  as  they  are  produced  at  a  low  pressure  helps 
maintain  this  r*j  <  situation,  and  rapid  expulsion  of  newly-created  electrons  keeps  ne  <C 
in  the  gas  while  n,-  <  n*  there.  If  the  cathode  is  hot  enough  to  emit  ionizing  radiation,  the 
gas  should  be  shielded  (or  shadowed)  from  it  in  order  to  avoid  photoionization. 

If,  as  at  higher  pressures,  n,-  builds  up  to  exceed  nj,  it  is  predicted  that  electron  density 
comparable  to  n<,  or  greater  can  build  up  in  the  source  gas,  leading  to  rapid  electron  heating 
by  two-stream  instability,  electron- impact  avalanche  of  the  remainder  of  the  gas,  and  a 
thermal  explosion  producing  a  denser  jet  of  plasma  with  much  larger  transverse  (thermal) 
energies  than  in  the  case  of  strictly  “electrostatic”  acceleration.  During  the  course  of  beam 
injection  at  pressures  on  the  order  of  1  Torr  and  above,  both  the  fast  tenuous  directed  jet 
and,  subsequently,  the  slower,  denser  jet  of  the  thermal  explosion  are  predicted,  and  are  seen 
in  a  ID  simulation.  After  beam  turnoff,  at  low  pressures  a  plasmoid  with  n,  ~  nj  is  seen  in 
simulations  to  separate  from  the  source  gas  and  propagate  downstream.  At  higher  pressures, 
simulation  shows  this  being  followed,  even  overtaken,  by  a  thermal-explosion  jet  in  which 
plasma  electrons  have  been  very  strongly  heated. 

The  interaction  of  a  propagating,  expanding  plasmoid  with  a  very  tenuous  ionized  back¬ 
ground  plasma  has  also  been  considered.  Electrostatic  effects  are  expected  to  dominate  over 
collisional  or  magnetic  effects.  Enhanced  electric  field  fluctuations  in  the  propagating  plas¬ 
moid  are  predicted  to  raise  its  transverse  electron  temperature  and  thus  contribute  to  its 
reduction  in  axially-directed  energy  flux  density  due  to  enhanced  radial  expansion. 
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Figure  Captions 


Figure  1 


Figure  2 


Figure  3 


Beam-electron,  plasma-ion,  and  plasma-electron  density  profiles  and  net-charge-density 
profile  at  time  steps  (a)  20,000,  (b)  22,000,  (c)  24,000,  and  (d)  26,000.  Number  densities 
are  shown  in  units  of  number  cm-3,  charge  densities  are  shown  in  units  of  C  m-3,  and 
distances  are  measured  from  the  anode  plane  in  cm. 

Beam-electron,  plasma-ion,  and  plasma-electron  density  profiles  and  net-charge-density 
profile  at  time  step  40,000  for  a  run  in  which  the  source  plasma  has  been  artificially 
removed. 

Momentum  distributions  for  (a)  beam  electrons  and  (b)  plasma  ions  at  20  cm  at  40,000 
time  steps  for  a  run  in  the  intermediate-pressure  regime. 


30 


100 


80 


60 


40 


nt  40000 

z  :  20.000 

A  :  10.000 

<u/c>  :  0.17932 

<y>  :  1.0524 

</?>  0.15676 


cm 
cm  - 


20 


nt 


jimitm 


z  : 

A  : 
<u/c>  : 
<7 >-/  : 

<P>  : 


20.000  cm  - 

10.000  cm 

3.45005E— 02 
6.00540E— 04 
3.44795E-02 


10  Papers  and  Presentations  Resulting  from  this  Research 

The  following  papers,  wholly  or  partly  derived  from  this  contract,  were  presented: 


1.  R.L.  Yao,  C.D.  Striffler  and  J.  Guillory,  “Intense  beam  propagation  and  collective  ion 
acceleration  in  a  Localized  plasma  source-vacuum  system,”  IEEE  Int  Conf.  on  Plasma 
Sci.,  June  1-3,  1987,  Crystal  City,  VA. 

2.  R.L.  Yao,  C.D.  Striffler  and  J.  Guillory,  “Collective  Ion  Acceleration,  Intense  Beam 
Propagation,  and  Plasmoid  Emission  in  a  Localized  Plasma  Source- Vacuum  System,” 
APS-DPP  Meeting,  Sand  Diego,  CA,  November  2-6,  1987  [Bull.  APS  32,  1719  (1987)]. 

3.  R.L.  Yao,  C.D.  Striffler  and  J.  Guillory,  “Numerical  Simulation  of  Intense  Beam  Prop¬ 
agation  and  Plasmoid  Formation  in  a  Localized  Plasma-Source  and  Vacuum  System,” 
IEEE  Int.  Conf.  on  Plasma  Science,  Seattle,  WA,  June  6-8,  1988;  also  APS-DPP 
Meeting,  Hollywood,  FL,  October  31-November  4,  1988  [Bull.  APS  33,  1955  (1988)]. 

4.  W.W.  Destler,  J.  Rodgers,  Z.  Segalov,  C.D.  Striffler,  R.L.  Yao,  J.  Guillory,  and  X. 
Zhang,  “Review  of  Intense  Electron  Beam  Transport  in  Gases,”  Proceedings,  Int.  Conf. 
BEAMS  ’88,  Karlsruhe,  W.  Germany,  July  4-8,  1988. 

5.  J.  Guillory,  R.L.  Yao,  C.D.  Striffler  and  M.  Reiser,  “University  of  Maryland  Plasmoid 
Program:  Formation  and  Stability  of  Partially- Neutralized  Plasma  Clumps  -  Plasmoid 
Formation  Concepts,”  DOE  Plasmoid  Workshop  at  SAIC,  McLean,  VA,  October  18-19, 
1988. 

6.  J.  Guillory,  “Interaction  of  a  Partially  Neutralized  Hot-Electron  Plasmoid  with  Low- 
Density  Background  Plasma  -  Preliminary  Theory,”  APS-DPP  Meeting,  Hollywood, 
FL,  November  4,  1988  [Bull.  APS  33,  2118  (1988)]. 

7.  J.  Guillory  and  R.L.  Yao,  “Theory  of  Ionization  of  Gas  by  a  Pulsed  Reflexing  Electron 
Beam  -  Loss  of  Stability,”  APS-DPP  Meeting,  Anaheim,  CA,  November  13-17,  1989 
[Bull.  APS  34,  2103  (1989)]. 


38 


8.  R.L.  Yao  and  C.D.  Striffler,  “Numerical  Simulation  of  Collective  Ion  Acceleration  in  an 
Intense  Electron  Beam  -  Localized  Gas  Cloud  System,”  J.  Appl.  Phys.  67,  1650-1658 
(1990).  [Copy  attached.] 

9.  J.  Guillory,  R.  Yao,  C.D.  Striffler,  and  M.  Reiser,  manuscript  in  final  stages  of  prepar- 
tion  for  submission  to  a  refereed  journal,  “Formation  and  Stability  of  Partially-Neutralized 
Plasma  Clumps.” 


39 


j 

Numerical  simulation  of  collective  ion  acceleration  in  an  intense  electron 
beam-localized  gas  cloud  system 

R.  L.  Yao  and  C.  D.  Striffler 

Laboratory  for  Plasma  Research  and  Department  of  Electrical  Engineering,  University  of  Maryland, 

College  Park,  Maryland  20742 

(Received  5  September  1989;  accepted  for  publication  7  November  1989) 

In  experiments  in  which  an  intense  relativistic  electron  beam  is  injected  into  an  evacuated  drift 
tube  with  a  localized  gas  cloud  located  near  the  anode,  ions  with  energies  several  times  the 
electron  beam  energy  have  been  observed.  These  experiments  have  been  simulated  using  a 
particle-in-cell  code  which  realistically  models  ionization  of  the  gas.  It  was  found  that  when 
the  injected  electron  beam  current  exceeds  the  space-charge  limiting  current,  ions  are 
accelerated  to  energies  several  times  the  electron  beam  energy  by  coherent  motion  of  the  ions 
and  the  intense  virtual  cathode  electric  fields.  The  dependence  of  the  peak  ion  energy  on  the 
system  parameters  as  observed  in  the  simulations  is  also  discussed.  For  the  parameter  regimes 
investigated  with  beam  energies  up  to  3  MV,  beam  currents  up  to  35  kA,  gas  pressures  up  to 
600  mTorr,  and  gas  cloud  widths  up  to  6  cm,  peak  ion  energies  of  5-6  times  the  electron  beam 
energy  have  been  observed. 


I.  INTRODUCTION 

Because  the  electric  fields  produced  by  an  intense  rela¬ 
tivistic  electron  beam  (IREB)  are  much  larger  than  those 
that  can  be  attained  in  conventional  accelerators,  consider¬ 
able  work  has  been  devoted  to  using  IREBs  to  accelerate  ions 
to  high  energies.  In  1970,  Graybill  and  Uglum 1  reported  that 
ions  with  peak  energies  of  four  times  the  electron  beam  ener¬ 
gy  were  produced  when  an  IREB  was  injected  into  a  gas- 
filled  drift  tube.  In  1973,  Luce,  Sahlin,  and  Crites2  reported 
that  ions  with  energies  several  times  the  electron  beam  ener¬ 
gy  were  produced  when  an  IREB  was  injected  into  an  evacu¬ 
ated  drift  tube  using  a  diode  with  a  dielectric  insert  in  the 
anode.  Since  these  results  were  published,  collective  ion  ac¬ 
celeration  has  been  reported  for  many  other  systems,  includ¬ 
ing  those  in  which  the  ions  are  produced  by  beam  ionization 
of  dielectric  drift  tube  walls1  or  by  laser  ionization  of  a  solid 
target.4 

Recently,  Destler,  Floyd,  and  Reiser’  reported  that  in 
experiments  performed  at  the  University  of  Maryland  col¬ 
lective  ion  acceleration  was  observed  when  an  IREB  was 
injected  into  an  evacuated  drift  tube  after  gas  had  been 
puffed  into  the  drift  tube  through  a  nozzle  located  near  the 
anode.  A  diagram  of  the  University  of  Maryland  experiment 
is  shown  in  Fig.  1 .  The  system  consists  of  a  foilless  electron 
beam  diode,  a  grounded  cylindrical  drift  tube,  a  nozzle  locat¬ 
ed  near  the  anode  through  which  gas  is  puffed  into  the  drift 
tube,  and  magnetic  field  coils,  which  produce  a  guide  field 
for  the  electron  beam.  In  a  typical  experiment,  the  drift  tube 
is  first  pumped  down  to  vacuum  (5x  10  5  Torr).  The  puff 
value  connected  to  the  nozzle  is  then  fired,  injecting  neutral 
gas  into  the  drift  tube  through  the  nozzle  and  producing  a 
cloud  of  gas  extending  no  more  than  a  few  cm  from  the  an¬ 
ode.  The  diode  is  then  fired,  injecting  the  beam  into  the  drift 
tube  through  the  cloud  of  gas.  Current  collectors  placed  at 
the  downstream  end  of  the  drift  tube  measure  the  ion  current 
and  stacked-foil  neutron  activation  analysis  is  used  to  deter¬ 
mine  the  ion  energy  and  flux. 

In  the  experiments,  it  was  found  that  when  the  injected 


current  70  exceeded  the  space-charge  limiting  current  IL  for 
the  system,  ions  were  accelerated  to  energies  several  times 
the  electron  beam  energy  V0.  For  example,  using  a  1.5-MeV, 
35-kA  electron  beam,  protons  were  accelerated  to  energies 
of  greater  than  8  MeV  (5.3  //nucleon )  and  xenon  ions 
were  accelerated  to  energies  of  greater  than  638  MeV  (4.9 
//nucleon  )f  It  was  also  found  that  the  maximum  ion  ener¬ 
gy  was  essentially  independent  of  the  applied  magnetic 
field.7 

The  space-charge  limiting  current  for  a  system  is  the 
maximum  current  that  can  be  propagated  through  the  sys¬ 
tem  before  the  potential  depression  produced  by  the  space- 
charge  exceeds  the  beam  energy.  When  the  injected  current 
exceeds  the  space-charge  limiting  current,  a  virtual  cathode 
forms  where  the  potential  equals  the  beam  energy.  The  frac¬ 
tion  of  the  injected  current  which  exceeds  the  space-charge 
limiting  current  is  reflected  by  the  virtual  cathode  back 
through  the  anode  or  deflected  to  the  drift  tube  wall.  For  a 
solid  beam  of  radius  Rb  in  a  cylindrical  drift  tube  of  radius 
Rw  with  a  large  focusing  magnetic  field,  the  space-charge 
limiting  current  is  given  by  the  formula" 


h 


17  000(yy3-  1)1/2 
1  +21n  (Rw/Rb) 


—  (Amps), 


(1) 


where  y0  is  the  Lorentz  factor  for  the  injected  electrons  and 
where  it  is  assumed  that  Rb  4RW- 

Many  explanations  have  been  offered  for  collective  ac¬ 
celeration  in  both  gas-filled  systems  and  evacuated  systems. 
Among  the  first  models  to  be  proposed  for  collective  ion 
acceleration  in  a  gas-filled  drift  tube  was  that  of  Olson.9  Ac¬ 
cording  to  this  model,  which  is  based  on  work  by  Poukey  and 
Rostoker, 10  collecti  ve  ion  acceleration  occurs  only  when  the 
current  exceeds  the  space-charge  limiting  current.  Soon 
after  the  electron  beam  is  injected,  the  beam  front  stops  near 
the  anode,  forming  a  stationary  potential  well  with  a  depth  of 
two  to  three  times  the  electron  beam  energy.  As  the  space 
charge  near  the  anode  is  neutralized  by  ionization  of  the  gas, 
the  beam  front  and  the  potential  well  suddenly  start  to  prop- 
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FIG.  1 .  Schematic  of  experimental  setup  for  University  of  Maryland  collec¬ 
tive  acceleration  experiments. 

agate  downstream.  Ions  are  accelerated  by  the  deep,  station¬ 
ary  potential  well  and  by  the  moving  potential  well.  This 
theory  accounts  for  many  of  the  experimental  observations 
of  collective  ion  acceleration  in  gas-filled  systems  and  is  con¬ 
sistent  with  some  2-D  particle  simulations  of  IREB  injection 
into  a  neutral  gas." 

The  energies  that  are  typically  observed  in  evacuated 
systems  such  as  the  Luce  diode  configuration  are  much  too 
high  to  be  explained  by  the  theory  of  Olson.  Furthermore, 
Adler,  Nation,  and  Serlin12  have  also  reported  detailed  ob¬ 
servations  of  collective  ion  acceleration  in  a  Luce  diode  con¬ 
figuration  which  indicate  that  ion  acceleration  was  not  pro¬ 
duced  by  a  deep  potential  well  or  by  a  moving  potential  well 
at  the  beam  front.  They  suggested  instead  that  ions  were 
accelerated  by  electrostatic  waves  which  are  driven  by  the 
two-stream  instability.  Destler  etal.4  suggested  that  in  evac¬ 
uated  systems  collective  acceleration  occurs  by  a  bootstrap 
process.  The  early  part  of  the  beam  pulse  forms  a  dense  plas¬ 
ma  near  the  anode.  If  the  beam  current  exceeds  the  space- 
charge  limiting  current  and  the  plasma  density  is  sufficiently 
high,  a  virtual  cathode  forms  at  the  downstream  edge  of  the 
plasma.  The  intense  virtual  cathode  electric  fields  draw  ions 
out  of  the  plasma,  the  ions  neutralize  the  electron  beam 
space  charge,  and  the  virtual  cathode  moves  further  down¬ 
stream.  As  this  process  is  repeated,  the  beam  electrons  trans¬ 
fer  momentum  to  the  ions  until  the  ions  attain  the  same  ve¬ 
locity  as  the  beam-front  electrons.  A  similar  model  was 
proposed  earlier  by  Ryutov  and  Stupakov11  for  collective 
ion  acceleration  produced  by  injection  of  an  IREB  through  a 
dense  plasma  layer.  Mako  and  Tajima14  modified  the  model 
using  data  obtained  from  experiments  with  a  gas- filled  drift 
tube,  and  were  able  to  recover  the  experimentally  measured 
peak  ion  energy  (ten  times  the  electron  beam  energy)  and 
ion  distribution.  They  also  used  a  particle  code  to  simulate 
injection  of  an  IREB  into  an  evacuated  drift  tube  with  a 
dense  plasma  layer  located  next  to  the  injection  plane.  The 
simulations  showed  that  the  peak  ion  energies  were  limited 
by  phase  instability  of  the  ions  and  the  accelerating  electric 
fields.  However,  using  an  artificially  low  ion  mass,  they  also 
found  that  the  peak  ion  energies  never  exceeded  three  times 
the  electron  beam  energy  and  that,  when  the  ion  mass  was 
doubled,  the  ion  velocity  decreased  inversely  as  the  square 
root  of  the  mass.  Using  a  2-1/2-D  electromagnetic  particle 
code  with  a  realistic  ion  mass,  Faehl  and  Peter15  and  Peter  et 
o/.16  simulated  the  same  system  and  found  that  the  virtual 
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cathode  formed  at  the  downstream  edge  of  the  plasma  layer 
soon  after  the  beam  was  injected  and  that  no  subsequent 
motion  of  the  virtual  cathode  occurred.  They  also  pointed 
out  that  no  such  motion  should  be  expected  because  the 
maximum  ion  current  that  can  be  extracted  from  the  plasma 
is  never  large  enough  to  neutralize  the  electron  beam  space 
charge.  They  suggested  that  collective  acceleration  is  the  re¬ 
sult  of  ambipolar  expansion  of  the  plasma  rather  than  mo¬ 
tion  of  the  virtual  cathode. 

In  simulations  of  injection  of  an  IREB  from  a  plasma 
into  free  space,  Chang  and  Reiser17  found  that  as  ions  were 
extracted  from  the  surface  of  the  plasma  by  the  virtual  cath¬ 
ode  electric  fields,  the  location  of  the  virtual  cathode  mini¬ 
mum  shifted  downstream  from  its  initial  location.  They  also 
observed  ions  which  had  been  accelerated  to  energies  ap¬ 
proximately  twice  the  electron  beam  energy  and  proposed 
that  the  coherent  motion  of  the  virtual  cathode  and  the  ions 
as  the  virtual  cathode  shifts  downstream  to  its  final  position 
was  responsible  for  accelerating  the  ions. 

In  this  paper  we  present  a  detailed  picture  of  the  accel¬ 
eration  mechanism  in  the  University  of  Maryland  experi¬ 
ments,  in  which  an  intense  relativistic  electron  beam  is  in¬ 
jected  into  an  evacuated  drift  tube  with  a  localized  gas  cloud 
located  near  the  anode.  We  have  obtained  this  picture  using 
particle-in-cell  (PIC)  code  simulations  which  realistically 
model  ionization  of  the  neutral  gas. 

Our  simulations  differ  from  previous  simulations  of  col¬ 
lective  acceleration  in  the  following  respects.  In  our  simula¬ 
tions,  an  IREB  is  injected  into  an  evacuated  drift  tube  with  a 
well-localized  gas  cloud  located  next  to  the  anode,  as  op¬ 
posed  to  a  gas-filled  drift  tube  or  an  evacuated  drift  tube  with 
a  plasma  layer  located  next  to  the  anode.  Ions  are  produced 
in  the  drift  tube  by  realistically  modeling  ionization  of  the 
gas  by  impact  ionization,  as  opposed  to  assuming  an  initial 
distribution  of  ions  or  injecting  ions.  Realistic  ion  masses  are 
also  used.  The  finite  radius  of  the  beam  is  also  taken  into 
account,  so  that  the  simulations  are  more  realistic  than 
strictly  1-D  simulations.  Finally,  the  IREB  is  injected  into  a 
grounded  drift  tube  rather  than  free  space  so  that  inductive 
effects  can  confidently  be  neglected.  (We  also  note  that  our 
code  is  simple  enough  to  permit  a  detailed  investigation  of 
the  acceleration  mechanism.) 

In  our  simulations,  we  find  that  when  the  injected  elec¬ 
tron  current  exceeds  the  space-charge  limiting  current,  a  vir¬ 
tual  cathode  forms  near  the  anode  plane  soon  after  the  beam 
is  injected.  Neutralization  of  the  electron  beam  space  charge 
by  impact  ionization  of  the  neutral  gas  then  causes  the  vir¬ 
tual  cathode  to  move  downstream  to  the  end  of  the  gas  re¬ 
gion,  where  the  virtual  cathode  stops.  As  the  virtual  cathode 
moves  through  the  gas,  ions  are  accelerated  to  energies  sev¬ 
eral  times  the  electron  beam  energy  by  the  coherent  motion 
of  the  ions  and  the  intense  virtual  cathode  electric  fields. 
(This  picture  is  similar  to  the  mechanism  observed  by  Chang 
and  Reiser. ) 

Before  proceeding,  we  would  like  to  point  out  that  the 
general  features  of  the  above  picture  are  similar  to  results 
which  have  been  obtained  earlier  at  the  University  of  Mary¬ 
land  . 1 8  The  present  work  has  also  been  extended  to  simulate 
enhanced  collective  ion  acceleration  in  a  laser-controlled 
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us,  drift  tube  length,  neutral  gas  pressure,  and  gas  cloud 
width.  We  present  a  picture  of  the  collective  acceleration 
mechanism  that  we  have  observed  in  all  of  the  simulations  in 
which  collective  acceleration  occurred.  To  do  so,  we  will  use 
the  results  of  a  single  simulation.  (We  note  that  the  param¬ 
eter  values  differ  from  typical  parameter  values  in  the  Uni¬ 
versity  of  Maryland  experiments.  However,  simulations  per¬ 
formed  with  the  experimental  parameter  values  show  that 
the  acceleration  mechanism  is  the  same  as  described  here. ) 

In  the  simulation,  a  15-kA,  500-keV,  0.5-cm-radius 
electron  beam  is  injected  into  a  4-cm-radius,  40-cm-long 
drift  tube  with  a  1-cm-wide,  450-mTorr  hydrogen  cloud  lo¬ 
cated  next  to  the  anode  plane.  The  beam  current  and  beam 
energy  are  assumed  to  rise  instantaneously.  According  to 
Eq.  ( 1 ),  the  space-charge  limiting  current  for  the  system  is 
approximately  1.4  kA.  The  time  step  is  approximately  0.75 
ps  and  the  grid  cell  width  is  approximately  0.2  mm. 

Figure  3  shows  the  time-integrated  proton  energy  spec¬ 
trum  7.5  cm  from  the  anode  plane  1 1  ns  after  the  start  of  the 
run.  The  energy  spectrum  can  be  roughly  divided  into  two 
parts — the  bulk  of  the  spectrum,  which  is  centered  about  0.6 
MeV,  and  a  tail  that  extends  out  to  1.5  MeV  (three  times  the 
electron  beam  energy ) .  We  ignore  the  ions  in  the  bulk  of  the 
spectrum  and  focus  on  the  ions  in  the  high-energy  tail.  Note 
that  since  the  potential  depression  at  the  center  of  the  drift 
tube  is  approximately  —  500  kV,  the  ion  spectrum  that 
would  be  observed  by  a  grounded  detector  would  be  the 
same  as  the  spectrum  shown  in  Fig.  3  shifted  to  the  left  by  0.5 
MeV.  In  the  following,  we  present  a  detailed  picture  showing 
how  the  ions  in  the  high-energy  tail  are  accelerated  to  ener¬ 
gies  exceeding  the  electron  beam  energy. 

In  Fig.  4,  the  potential  <S>  on  axis  near  the  anode  plane  is 
plotted  as  a  function  of  distance  z  from  the  anode  plane  at 
several  different  time  steps.  About  40  ps  after  the  start  of  the 
run,  a  virtual  cathode  forms  about  2  mm  from  the  anode 
plane,  and  the  current  measured  downstream  from  the  vir¬ 
tual  cathode  falls  to  a  value  approximately  equal  to  the 
space-charge  limiting  current.  After  formation  of  the  virtual 
cathode,  the  potential  remains  essentially  the  same  for  sever¬ 
al  thousand  time  steps  (the  position  and  amplitude  of  the 
virtual  cathode  oscillate  in  time  at  the  local  plasma  frequen- 
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FIG.  3.  Time-integrated  proton  energy  spectrum  at  z  =  7.5  cm  after  1 1  ns. 
The  run  parameters  are  /„  =  15  kA,  V„  =  500  keV.  Rh  =  0.5  cm,  Rn  =  4 
cm,  d  =  40  cm.  p„  450  mTorr,  and  z„  =  1  cm. 
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FIG.  4.  Potential  on  axis  4>  near  anode  plane  vs  distance  from  anode  plane  z 
for  time  steps  4600, 4800,  5000,  5200,  and  5600  during  acceleration  phase. 
Note  that  1000  time  steps  equals  about  0.75  ns.  See  Fig.  3  for  system  param¬ 
eters. 


cy,  but  the  time-averaged  position  and  amplitude  are  ap¬ 
proximately  the  same).  The  potential  for  time  step  4600, 
which  is  shown  in  Fig.  4,  is  typical  for  this  stage.  During  this 
stage,  the  beam  electrons  ionize  the  neutral  gas,  producing  a 
plasma  in  the  vicinity  of  the  virtual  cathode.  The  intense 
electric  fields  near  the  virtual  cathode  quickly  expel  the  plas¬ 
ma  electrons  from  the  plasma  region.  The  plasma  ions  are 
also  accelerated  by  these  electric  fields  but  leave  the  plasma 
region  at  a  much  slower  rate  than  the  plasma  electrons.  They 
also  tend  to  be  accelerated  downstream,  whereas  the  plasma 
electrons  tend  to  be  accelerated  upstream,  where  they  exit 
from  the  drift  tube  through  the  anode  plane.  As  plasma  ions 
drift  through  the  gas,  they  also  contribute  to  the  ionization. 
As  the  ion  density  in  the  plasma  region  increases,  the  virtual 
cathode  moves  downstream  to  the  end  of  the  plasma  region 
leaving  behind  a  positive  potential  region  (time  steps  4600- 
5000) .  Ions  lying  to  the  left  of  the  peak  of  the  positive  poten¬ 
tial  are  accelerated  upstream  and  exit  the  drift  tube  through 
the  anode  plane.  Ions  lying  to  the  right  of  the  peak  of  the 
positive  potential  are  accelerated  downstream.  As  the  ions 
move  downstream,  the  virtual  cathode  also  moves  down¬ 
stream  with  the  positive  potential  region  widening  and  ex¬ 
tending  downstream  behind  it  (time  steps  5000-5600).  The 
virtual  cathode  ceases  to  move  when  it  reaches  the  end  of  the 
gas  region  at  z  =  1cm  (time  step  5600).  We  will  see  that  ions 
are  accelerated  to  energies  of  nearly  1.5  MeV  during  the  time 
that  the  virtual  cathode  and  the  positive  potential  region 
move  downstream. 

The  formation  of  a  positive  potential  region  has  been 
observed  experimentally  and  in  other  simulations.  I‘‘  The 
reason  the  potential  in  the  plasma  region  becomes  positive 
can  be  understood  as  follows.  Because  the  electron  impact 
ionization  cross  section  peaks  at  low  energies,  the  majority  of 
ions  are  formed  near  the  virtual  cathode.  The  buildup  of  ions 
in  the  vicinity  of  the  virtual  cathode  causes  the  beam  elec¬ 
trons  to  spend  more  time  there,  and  as  a  result  the  electron 
density  downstream  from  the  virtual  cathode  decreases. 
Hence,  when  the  ion  density  is  sufficiently  high,  the  position 
of  the  virtual  cathode  can  no  longer  be  maintained  and  the 
virtual  cathode  begins  to  move  downstream.  ( In  fact,  our 
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FIG.  5.  Phase-space  trajectory  of  collectively  accelerated  proton.  The  boxes 
plotted  along  the  trajectory  indicate  the  location  of  the  ion  at  equally  spaced 
time  steps  during  the  acceleration  phase.  The  diamonds  plotted  along  the 
base  of  the  trajectory  indicate  the  location  of  the  peak  electric  field  at  the 
same  time  steps.  See  Fig.  3  for  system  parameters. 


simulations  show  that  immediately  preceding  the  onset  of 
the  acceleration  process,  the  electron  current  measured  im¬ 
mediately  downstream  from  the  plasma  region  falls  abrupt¬ 
ly.  )  As  the  virtual  cathode  moves  through  the  plasma  re¬ 
gion,  the  electrons  behind  the  virtual  cathode  travel  through 
the  plasma  with  essentially  the  same  speed  with  which  they 
are  injected.  The  increase  in  the  speed  of  the  electrons 
through  the  plasma  region  produces  a  sudden  drop  in  the 
electron  density.  Hence,  the  virtual  cathode  leaves  a  net  posi¬ 
tive  space  charge  in  its  wake  and  the  potential  becomes  posi- 
^  tive  in  the  plasma  region. 

In  Fig.  5,  the  phase-space  trajectory  of  a  collectively 
accelerated  proton  is  shown  during  the  collective  accelera¬ 
tion  process.  The  ion  is  created  with  zero  velocity  at  a  dis¬ 
tance  of  about  1.3  mm  from  the  anode  plane  about  2.6  ns 
after  the  start  of  the  run  ( time  step  3500) .  It  is  rapidly  accel¬ 
erated  by  the  electric  field  near  the  virtual  cathode  to  a  veloc¬ 
ity  of  about  0.008c  (corresponding  to  an  energy  of  30  keV) 
and  begins  to  drift  downstream.  When  it  has  drifted  about  2 
mm  downstream,  the  potential  begins  to  go  positive  (at 
around  time  step  4700)  and  the  virtual  cathode  begins  to 
move  downstream.  As  the  virtual  cathode  moves  down¬ 
stream  and  the  positive  region  extends  downstream  behind  it 
(time  steps  4800-5600),  the  proton  is  rapidly  accelerated  to 
a  final  velocity  of  about  0.054c  (corresponding  to  an  energy 
of  1.4  MeV). 

The  boxes  plotted  along  the  proton  phase-space  trajec¬ 
tory  indicate  points  along  the  trajectory  which  are  separated 
by  equal  time  intervals  of  200  time  steps  or  about  1 50  ps. 
Each  box  is  labeled  with  the  corresponding  time  step  starting 
with  time  step  4800.  Below  the  phase-space  trajectory,  we 
have  also  plotted  the  location  of  the  peak  electric  field  at  the 
same  time  steps.  The  peak  electric  field  is  always  located  in 
the  region  where  the  potential  drops  sharply  from  the  edge  of 
the  positive  potential  region  to  the  virtual  cathode.  As  the 
virtual  cathode  moves  downstream  and  the  positive  poten¬ 
tial  region  extends  downstream  behind  the  virtual  cathode 
the  peak  electric  field  also  moves  downstream.  As  can  be 
seen  from  the  figure,  the  ion  follows  closely  behind  the  peak 


electric  field  for  time  steps  4800  to  5600.  It  is  this  coherent 
motion  of  the  ion  and  the  peak  electric  field  through  the 
plasma  region  which  enables  the  ion  to  be  accelerated  to 
energies  several  times  the  energy  it  would  gain  if  the  poten¬ 
tial  were  stationary. 

Note  that  the  trajectory  of  the  proton  in  Fig.  5  is  similar 
to  the  trajectory  of  an  ion  being  accelerated  by  a  constant 
electric  field  vz(z)  =  y[v^  +  2 (eE2/mp)z,  where  vM  is  the 
initial  ion  velocity.  Using  values  obtained  from  the  figure  in 
the  above  formula,  we  find  that  the  electric  field  Ez  is  about 
200  MV/m,  which  is  near  the  peak  virtual  cathode  electric 
field  of  250  MV/m. 

The  reason  for  the  coherent  motion  of  the  ion  and  the 
peak  electric  field  during  the  motion  of  the  potential  front 
can  be  understood  by  considering  the  free  expansion  of  a  slab 
of  ions,  which  initially  has  uniform  density.  As  the  slab  ex¬ 
pands,  an  ion  at  the  surface  of  the  slab  always  feels  the  same 
electric  field,  that  is,  the  ion  and  the  electric  field  travel  at  the 
same  velocity.  We  believe  that  the  same  effect  occurs  during 
motion  of  the  potential  front.  In  the  simulations,  the  ions  at 
the  front  of  the  accelerated  ion  beam  are  accelerated  not  only 
by  the  space  charge  produced  by  the  ions  following  behind 
but  also  by  the  electric  field  produced  by  the  virtual  cathode 
which  always  moves  in  front  of  the  ion  beam — the  ions  are 
pushed  ahead  by  the  ions  which  follow  and  simultaneously 
pulled  ahead  by  the  virtual  cathode. 

The  reason  the  potential  front  stops  moving  at  the  end  of 
the  gas  region  can  also  be  understood  by  considering  the  free 
expansion  of  a  slab  of  ions.  As  the  slab  expands,  the  ion 
density  and  the  potential  decrease  uniformly.  In  the  gas  re¬ 
gion,  ionization  of  the  gas,  particularly  by  ions,  can  supply 
ions  at  a  rate  sufficiently  high  to  prevent  total  collapse  of  the 
potential  and  reformation  of  a  virtual  cathode  near  the  an¬ 
ode.  Beyond  the  gas  region,  no  ions  are  produced  and  the  ion 
density  is  too  low  to  permit  the  virtual  cathode  to  move  for¬ 
ward.  (As  Faehl  and  Peter15  observed,  the  maximum  ion 
current  which  can  be  extracted  from  the  plasma  is  space- 
charge  limited  and  hence  ions  can  never  be  extracted  at  a  rate 
sufficient  to  significantly  neutralize  the  electron  beam  space 
charge  downstream.) 

Although  all  of  the  collectively  accelerated  ions  are  ac¬ 
celerated  by  the  moving  potential  front,  some  ions  continue 
to  accelerate,  although  at  a  much  lower  rate,  after  they  have 
left  the  gas  region.  This  additional  acceleration  is  probably 
produced  by  the  interaction  of  the  collectively  accelerated 
ions  among  themselves — interactions  between  ions  that 
have  been  collectively  accelerated  in  a  bunch  will  cause  ions 
with  velocities  greater  than  the  mean  velocity  to  be  acceler¬ 
ated  and  ions  with  velocities  less  than  the  mean  velocity  to  be 
decelerated. 

Several  bursts  of  accelerated  ions,  which  follow  one  an¬ 
other  at  short  intervals,  are  often  observed  in  simulations. 
The  multiple  bursts  are  produced  when  the  positive  potential 
region  collapses  after  the  ions  have  been  accelerated  out  of 
the  gas  region,  the  virtual  cathode  retreats  to  a  location  in¬ 
side  the  gas  region  and  the  acceleration  process  repeats  itself. 
If  the  pressure  is  increased,  with  all  other  system  parameters 
kept  constant,  the  multiple  bursts  are  suppressed,  probably 
because  the  rate  of  ion  production  in  the  gas  region  is  suffi- 
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cient  to  prevent  collapse  of  the  positive  potential  region. 
Such  multiple  bursts  have  been  observed  in  some  experi¬ 
ments  in  evacuated  systems. 12 

IV.  RESULTS  OF  PARAMETER  STUDIES 
A.  Dependence  on  beam  energy  and  current 

In  Fig.  6,  the  peak  proton  energy  EmdX  measured  at 
z  =  7.5  cm  in  the  first  10  ns  of  the  run  is  plotted  as  a  function 
of  the  injected  beam  current  70  for  beam  energies  V0  of  0.5, 
1.0,  1.5,  and  3.0  MeV.  In  each  case,  the  gas  pressure  is  600 
mTorr  and  the  gas  cloud  width  is  1  cm.  Generally,  £max 
increases  with  beam  current  and  beam  energy.  In  Fig.  7,  the 
same  data  are  plotted  with  normalized  to  V0  and  70 
normalized  to  the  space-charge  limiting  current  IL ,  given  by 
Eq.  (1).  Notice  that  all  the  plotted  points  lie  on  approxi¬ 
mately  the  same  curve.  For  all  values  of  V0  which  were  used, 
no  collectively  accelerated  protons  are  observed  for 
7</7t  <  1.  As  IJIL  increases  above  unity,  increases 
sharply  and  then  continues  to  increase  but  at  a  much  slower 
rate.  (Note  that  for  fixed  7„,  the  ratio  Io/IL  decreases  with 
increasing  V0.  Hence  for  the  same  current  70,  the  ratio 
E„../Vn  with  V0  =  0.5  MeV  is  much  larger  than  the  ratio 
for  V0  =  3  MeV.) 

The  observed  dependence  of  the  peak  proton  energy  on 
beam  energy  and  current  can  be  understood  q»'.ke  simply  by 
recalling  that  the  strength  of  the  peak  electric  field  increases 
with  V0  and  IJ1L  ■  In  Fig.  8(a),  the  phase-space  trajectories 
of  collectively  accelerated  protons  are  shown  for  three  differ¬ 
ent  values  of  the  injected  current — 10,  20,  and  30  kA — and 
the  same  beam  energy  of  0.5  MeV.  In  each  case,  the  gas 
pressure  is  600  mTorr  and  the  gas  cloud  width  is  1  cm.  Re¬ 
calling  that  the  accelerating  electric  field  is  approximately 
proportional  to  the  product  of  the  velocity  u  and  the  slope 
du/dz  of  the  trajectory,  it  can  be  seen  that  as  the  current 
increases  so  does  the  accelerating  electric  field.  In  Fig.  8(b), 
the  phase-space  trajectories  of  collectively  accelerated  pro¬ 
tons  are  shown  for  three  different  values  of  the  beam  ener¬ 
gy — 0.5,  1.0,  and  1.5  MeV — and  the  same  beam  current  of 
30  kA.  It  can  be  seen  that  as  the  beam  energy  increases  so 
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FIG.  6.  Peak  proton  energy  Emm,  measured  at  z  =  7.5  cm  in  the  first  10  ns  vs 
injected  current  /„  for  p„  =  600  mTorr,  z„  =  1  cm,  and  F„  =  0.5,  1.0,  1.5, 
and  3.0  MeV.  See  Fig.  3  for  system  parameters. 


FIG.  7.  Peak  proton  energy  EmM ,  normalized  to  beam  energy  Va  vs  injected 
current  /0  normalized  to  space-charge  limiting  current  IL .  See  Figs.  3  and  6 
for  system  parameters. 


does  the  accelerating  electric  field.  (This  dependence  is  not 
so  clear  when  the  trajectories  for  0.5  and  1.0  MeV  are  com¬ 
pared,  which  illustrates  the  complex  nature  of  the  accelera¬ 
tion  process,  that  is,  the  final  ion  energy  depends  not  only  on 
the  strength  of  the  electric  field  but  also  on  the  location  and 
velocity  of  the  ion  when  the  virtual  cathode  begins  to  move 
downstream  and  on  the  velocity  with  which  the  virtual  cath¬ 
ode  moves. ) 
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FIG.  8.  Phase-space  trajectories  of  collectively  accelerated  protons  for  (a) 
V„  =  0.5  MeV  and  /„  =  10.  20.  and  30  kA,  and  for  (b)  /„  =  30  kA  and 
Vu  =  0.5,  1.0,  and  1.5  MeV.  See  Fig.  3  for  other  system  parameters. 
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(C)  p0  (Torr) 

FIG.  9.  Peak  proton  energy  measured  at  z  —  7.5  cm  after  20  ns  vs  gas 
pressure p„  for  z„  =  1  cm  and  K„  =  (a)  0.5,  (b)  1.0,  and  (c)  1.5  MeV.  See 
Fig.  3  for  other  system  parameters. 


B.  Dependence  on  gas  cloud  pressure 

In  Fig.  9,  the  peak  proton  energy  £max  measured  at 
z  =  7.5  cm  in  the  first  20  ns  of  the  beam  pulse  is  plotted  as  a 
function  of  gas  pressure  p0  for  injected  beam  currents  /„  of 
20,  25,  and  30  kA  and  for  beam  energies  V0  of  (a)  0.5,  (b) 
1 .0,  and  ( c )  1 . 5  Me V .  I  n  each  case,  the  gas  cloud  width  is  1 .0 
cm.  At  50  mTorr,  no  collectively  accelerated  protons  are 
observed  in  the  first  20  ns  of  the  beam  pulse.  Above  50 
mTorr,  Emax  tends  to  increase  with  increasing  pressure  and 
level  off  at  energies  several  times  the  electron  beam  energy. 
Generally,  the  higher  the  beam  energy,  the  higher  the  energy 
and  pressure  at  which  the  energy  levels  off. 


The  observed  dependence  of  the  peak  proton  energy  on 
the  pressure  can  be  explained  qualitatively  as  follows.  At  low 
pressures,  the  amount  of  ionization  at  the  end  of  20  ns  is  too 
small  to  neutralize  the  electron  beam  space  charge  and  the 
virtual  cathode  remains  at  its  initial  position;  no  collective 
acceleration  occurs  and  the  peak  proton  energies  which  are 
measured  are  less  than  the  electron  beam  energy.  As  the 
pressure  is  increased,  the  amount  of  ionization  which  is  pro¬ 
duced  during  the  beam  pulse  becomes  large  enough  to  neu¬ 
tralize  the  electron  beam  space  charge  and  permit  the  virtual 
cathode  to  move  downstream.  As  it  moves  downstream,  the 
virtual  cathode  electric  fields  accelerate  ions  to  energies  sev¬ 
eral  times  the  electron  beam  energy.  Once  the  acceleration 
process  has  begun,  the  final  energy  of  a  collectively  acceler¬ 
ated  ion  depends  only  on  the  magnitude  of  the  accelerating 
electric  field  and  the  distance  over  which  the  ion  travels  co¬ 
herently  with  the  electric  field.  Now,  the  magnitude  of  the 
electric  field  will  be  determined  mainly  by  the  electron  beam 
energy  and  current  and  the  density  of  ions  in  the  plasma 
region  at  the  start  of  the  acceleration  process  (because  the 
electric  field  at  the  surface  of  a  slab  of  ions  is  proportional  to 
the  ion  density  in  the  slab) .  But  the  ion  density  at  the  start  of 
the  acceleration  process  is  independent  of  the  pressure,  be¬ 
cause  it  is  proportional  to  the  electron  beam  density,  which  is 
determined  by  the  electron  beam  energy  and  current  only. 
The  only  effect  of  increasing  the  gas  pressure  will  be  to  de¬ 
crease  the  time  required  for  the  ion  density  to  build  up  before 
the  acceleration  process  begins  (which  is  in  fact  observed  in 
simulations).  Hence,  the  peak  ion  energy  will  depend  only 
weakly  on  the  neutral  gas  pressure,  once  the  gas  pressure  is 
sufficiently  high  to  cause  collective  acceleration. 

The  above  picture  is  somewhat  too  simplistic  because 
the  gas  pressui  also  affects  the  motion  of  the  virtual  cathode 
electric  fields  through  the  gas  cloud.  Generally,  as  the  pres¬ 
sure  increases,  the  virtual  cathode  moves  through  the  gas 
cloud  more  quickly  and  higher  peak  proton  energies  are  ob¬ 
served  because  the  virtual  cathode  electric  fields  are  able  to 
accelerate  the  protons  over  a  longer  distance. 

Although  the  dependence  of  collective  acceleration  on 
pressure  as  observed  in  our  simulations  is  consistent  with  the 
model  for  collective  acceleration  which  we  have  presented 
above,  it  disagrees  with  the  experimental  results  obtained  at 
the  University  of  Maryland.  In  particular,  Floyd  et  al.b  re¬ 
ported  that  when  a  1.5-MeV,  35-kA,  30-ns  beam  pulse  was 
injected  into  a  25-cm-diam  drift  tube  with  a  well-localized 
gas  cloud  located  near  the  anode,  collective  acceleration  was 
observed  only  in  a  narrow  pressure  window.  In  our  simula¬ 
tions,  we  observe  no  such  pressure  window.  In  particular,  we 
observe  collective  acceleration  at  all  pressures  up  to  1  Torr. 

The  disagreement  may  hr  /e  the  following  causes.  First, 
in  our  simulations  we  have  neglected  charge  exchange  be¬ 
tween  protons  and  H:  molecules.  Since  the  mean  free  path 
for  charge  exchange  for  protons  in  H2  is  approximately  1  cm 
at  100  mTorr  at  approximately  100  keV,9  at  pressures  of  a 
few  100  mTorr  charge  exchange  probably  prevents  accelera¬ 
tion  of  protons  to  energies  much  greater  than  a  few  100  keV. 
Second,  in  our  simulations  we  have  neglected  sources  of  ioni¬ 
zation  such  as  photoionization  and  prepulse  ionization 
which  could  produce  sufficient  ionization  in  the  gas  region 
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to  prevent  formation  of  a  virtual  cathode  near  the  injection 
plane.  In  simulations  of  injection  of  an  IREB  into  vacuum 
with  a  localized  plasma  located  next  to  the  injection 
plane,1516,21  it  was  found  that  for  plasma  densities  compara¬ 
ble  to  or  exceeding  the  beam  density,  the  virtual  cathode 
formed  near  the  downstream  edge  of  the  plasma.  Finally,  the 
disagreement  with  the  reported  experiments  may  be  the  re¬ 
sult  of  changes  in  the  electron  beam  diode  characteristics  as 
the  pressure  is  increased.  In  the  experiments,  the  peak  gas 
cloud  pressure  at  the  time  of  beam  injection  was  varied  by 
simply  varying  the  time  delay  between  firing  of  the  gas  puff 
valve  and  firing  of  the  diode.  It  was  also  found  that  when  the 
gas  puff  was  fired  too  early,  the  diode  shorted  out.4 


o  100  mTorr 
6  150  mTorr 
°  200  mTorr 


o  100  mTorr 
5  *  150  mTorr 

o  200  mTorr 


C.  Dependence  on  gas  cloud  width 

In  Fig.  10,  the  peak  proton  energy  Emii>  measured  at 
z  =  7.5  cm  in  the  first  20  ns  of  the  beam  pulse  is  plotted  as  a 
function  of  gas  cloud  width  z0  for  electron  beam  energies  V() 
of  (a)  0.5,  (b)  1.0,  and  (c)  1.5  MeV  and  several  reutral  gas 
pressures p0.  In  each  case,  the  injected  electron  beam  current 
is  30  kA. 

Figure  10  shows  that  £mas  tends  to  increase  with  gas 
cloud  width  and  then  level  off,  and  that  this  tendency  be¬ 
comes  more  pronounced  as  the  pressure  increases.  In  Fig. 
10(c),  for  example,  which  shows  data  obtained  for  V„  =  1.5 
MeV,  £max  increases  approximately  linearly  with  increasing 
gas  cloud  width  before  leveling  off  at  about  6.0  MeV  at 
z0  =  5.0  cm  for  p0  =  200  mTorr.  (When/;,,  =  100  mTorr,  no 
collectively  accelerated  protons  are  detected  in  the  first  20  ns 
of  the  beam  pulse. ) 

The  tendency  of  £max  to  increase  with  increasing  gas 
cloud  width  is  consistent  with  the  simple  picture  of  collective 
acceleration  presented  above.  Protons  are  accelerated  to 
higher  energies  as  the  gas  cloud  width  is  increased  because 
they  are  accelerated  over  longer  distances.  This  trend  is 
more  pronounced  at  higher  pressures  because,  as  noted  in 
the  last  section,  higher  pressures  enable  the  virtual  cathode 
to  move  through  the  gas  cloud  more  quickly.  Closer  exami¬ 
nation  of  the  simulation  results  shows  that  Emax  tends  to 
increase  as  the  gas  cloud  width  increases  because  the  protons 
are  indeed  accelerated  over  longer  distances.  Figure  11 
shows  the  phase-space  trajectories  of  two  collectively  accel¬ 
erated  ions  for  the  same  beam  energy  and  current  (30  kA) 
and  the  same  gas  cloud  pressure  (200  mTorr)  but  different 
gas  cloud  widths.  It  is  clear  from  the  figure  that  one  proton  is 
accelerated  to  a  higher  energy  because  it  is  accelerated  over  a 
longer  distance.  The  figure  also  illustrates  the  reason  for  the 
tendency  of  the  peak  proton  energy  to  level  off  as  the  gas 
cloud  width  is  increased— the  ions  tend  to  outrun  the  mov¬ 
ing  potential  front  as  it  moves  through  the  gas  cloud  (seen  in 
the  3 -cm  case). 

Results  of  simulations  performed  with  different  values 
of  the  injected  beam  current  /„  indicate  that  the  value  at 
which  £max  levels  off  also  tends  to  increase  with  beam  cur¬ 
rent,  although  the  correlation  is  relatively  weak.  In  Fig.  12, 
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FIG.  10.  Peak  proton  energy  E„,,x  measure  at  z  =  7.5  cm  after  20  ns  vs  gas 
cloud  width  2, ,  for/, |  =  30kA,p„  =  100,  150,  and  200  mTorr,  and  V„=  (a) 
0.5,  (b)  1.0,  and  (c)  1.5  MeV.  See  Fig.  3  for  other  system  parameters. 


FIG.  11.  Phase-space  trajectory  of  collectively  accelerated  proton  for 
V„  =  1.5  MeV,  /„  =  30kA./>„  =  200  mTorr,  and  2,,  =  2  and  3  cm.  See  Fig.  3 
for  other  system  parameters. 
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FIG.  12.  Peak  proton  energy  measured  at  z  =  7.5  cm  after  20  ns  vs  gas 
cloud  width  z„  for  V„  =  1.5  MeV.p,,  =  200  mTorr,  and  /„  =  20,  25,  and  30 
kA.  See  Fig.  3  for  other  system  parameters. 


for  example,  the  peak  proton  energy  measured  at  z  =  7.5  cm 
in  the  first  20  ns  of  the  beam  pulse  is  plotted  as  a  function  of 
gas  cloud  width  for  a  beam  energy  of  1 .5  MeV,  a  gas  pressure 
of  200  mTorr,  and  injected  beam  currents  of  20,  25  and  30 
kA. 

V.  SUMMARY 

In  order  to  investigate  collective  ion  acceleration  ob¬ 
served  in  experiments  in  which  an  intense  relativistic  elec¬ 
tron-beam  is  injected  into  an  evacuated  drift  tube  with  a 
localized  gas  cloud  located  near  the  anode,  we  have  written  a 
relativistic,  electrostatic  PIC  code  that  realistically  models 
ionization  of  the  gas. 

Using  the  code,  we  have  developed  a  simple  model  for 
the  collective  acceleration  mechanism  in  these  e  meriments. 
In  this  model,  ionization  of  the  gas  causes  the  virtual  cathode 
to  move  downstream  through  the  gas  cloud.  As  the  virtual 
cathode  moves  through  the  gas  cloud,  ions  are  accelerated  to 
energies  several  times  the  electron-beam  energy  by  coherent 
motion  of  the  ions  and  the  intense  virtual  cathode  electric 
fields. 

Using  the  code,  we  have  also  determined  the  depen¬ 
dence  of  the  peak  ion  energy  on  the  system  parameters,  in¬ 
cluding  the  beam  energy  V0,  beam  current  /„,  the  gas  cloud 
pressure p„,  and  the  gas  cloud  width  r„.  We  find  that  the  peak 
energy  tends  to  increase  with  increasing  values  of  beam 

energy  T()  and  beam  current  /„.  In  particular  we  find  that 
Ema%/Va  tends  to  increase  as  /,//,,  where  I,  is  the  space- 
charge  limiting  current  given  by  Eq.  ( 1 ),  and  that  the  func¬ 
tional  dependence  is  nearly  the  same  for  a  wide  range  of 
parameters.  We  also  find  that  as  the  pressure  is  increased, 
^max  jumps  from  a  value  approximately  equal  to  the  beam 
energy  to  a  value  several  times  the  beam  energy,  and  that 


generally,  the  higher  the  beam  energy,  the  higher  the  pres¬ 
sure  at  which  EmilK  reaches  the  final  energy.  We  also  find 
that,  generally,  Emax  tends  to  increase  with  gas  cloud  width 
z0  and  then  level  off  when  z0  is  a  few  cm  and  that  this  tenden¬ 
cy  is  more  pronounced  for  higher  pressures.  Our  results  indi¬ 
cate  that  the  mechanism  for  collective  acceleration  which  we 
have  observed  in  our  simulations  can  produce  protons  with 
energies  of  at  least  5-6  times  the  electron  beam  energy  when 
the  beam  and  the  gas  cloud  properties  are  optimized. 
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